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Abstract We review and discuss horizontal branch 
(HB) stars in a broad astrophysical context, includ- 
ing both variable and non-variable stars. A reassess- 
ment of the Oosterhoff dichotomy is presented, which 
provides unprecedented detail regarding its origin and 
systematics. We show that the Oosterhoff dichotomy 
and the distribution of globular clusters in the HB 
morphology-metallicity plane both exclude, with high 
statistical significance, the possibility that the Galac- 
tic halo may have formed from the accretion of dwarf 
galaxies resembling present-day Milky Way satellites 
such as Fornax, Sagittarius, and the LMC — an argu- 
ment which, due to its strong reliance on the ancient 
RR Lyrae stars, is essentially independent of the chem- 
ical evolution of these systems after the very earli- 
est epochs in the Galaxy's history. Convenient ana- 
lytical fits to isochrones in the HB type-[Fe/H] plane 
are also provided. In this sense, a rediscussion of the 
second-parameter problem is also presented, focusing 
on the cases of NGC 288/NGC 362, M13/M3, the ex- 
treme outer-halo globular clusters with predominantly 
red HBs, and the metal-rich globular clusters NGC 6388 
and NGC 6441. The recently revived possibility that 
the helium abundance may play an important role as a 
second parameter is also addressed, and possible con- 
straints on this scenario discussed. We critically discuss 
the possibility that the observed properties of HB stars 
in NGC 6388 and NGC 6441 might be accounted for 
if these clusters possess a relatively minor population 
of helium-enriched stars. A technique is proposed to 
estimate the HB types of extragalactic globular clus- 
ters on the basis of integrated far-UV photometry. The 
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importance of bright type II Cepheids as tracers of 
faint blue HB stars in distant systems is also empha- 
sized. The relationship between the absolute V mag- 
nitude of the HB at the RR Lyrae level and metallic- 
ity, as obtained on the basis of trigonometric parallax 
measurements for the star RR Lyr, is also revisited. 
Taking into due account the evolutionary status of RR 
Lyr, the derived relation implies a true distance mod- 
ulus to the LMC of (m-M)o = 18.44 ± 0.11. Tech- 
niques providing discrepant slopes and zero points for 
the My (RRL) — [Fe/H] relation are briefly discussed. 
We provide a convenient analytical fit to theoretical 
model predictions for the period change rates of RR 
Lyrae stars in globular clusters, and compare the model 
results with the available data. Finally, the conductive 
opacities used in evolutionary calculations of low-mass 
stars are also investigated. 

Keywords Galaxies: Local Group ■ Galaxy: forma- 
tion • Galaxy: globular cluster: general ■ Stars: evo- 
lution • Stars: Hertzsprung-Russell diagram • Stars: 
horizontal-branch • Stars: variables: other 



1 Introduction 

1.1 A Bit of History 

In her beautiful review of (hot) horizontal-branch (HB) 
stars, Moehler (2001) notes that Barnard (1900) was 
the first to detect the presence of (blue) horizontal- 
branch stars in globular clusters. The term horizontal 
branch appears to have been coined by ten Bruggen- 
cate (1927), to whom Moehler (2004) assigns the dis- 
covery of the horizontal branch — which he noticed when 
plotting the color-magnitude data obtained by Shapley 
(1915) in the latter's study of NGC 5272 (M3). Of 
course, with the development of nuclear astrophysics 
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and the establishment of modern stellar evolution the- 
ory still several years away, it was not until three 
decades later that Hoylc & Schwarzschild (1955) first 
correctly identified HB stars as the progeny of low-mass 
red giant branch (RGB) stars, burning helium in their 
center and hydrogen in a shell around the core. 

The first successful HB models were actually com- 
puted by Faulkner (1966), and Castellani & Renzini 
(1968) and Iben & Rood (1970) were the first to recog- 
nize that substantial mass loss on the RGB phase was 
needed to explain the observed colors of HB stars in 
globular clusters, with moreover a significant spread in 
mass loss amounts from star to star in any given globu- 
lar cluster being needed to explain their observed color 
ranges — blue HB stars losing, on average, more mass 
than red HB stars. The distribution of masses along 
the HB often resembles a normal or Gaussian distribu- 
tion (Rood & Crocker 1989; Dixon et al. 1996; Valcarce 
& Catelan 2008), and normal deviates are accordingly 
often adopted in the construction of "synthetic horizon- 
tal branches" (e.g., Rood 1973; Castellani & Tornambc 
1981; Caputo et al. 1987; Catelan 1993; Lee 1990; Lee, 
Demarque, & Zinn 1990; Cassisi et al. 2004). The pres- 
ence of mass distributions that resemble Gaussian de- 
viates strongly suggests the presence of stochastic mass 
loss processes on the RGB. However, deviations from 
a Gaussian shape are also not uncommon among glob- 
ular clusters, particularly in the cases of those having 
bimodal HBs and/or long blue tails with gaps (Cate- 
lan et al. 1998; Ferraro et al. 1998; Piotto et al. 1999; 
Momany et al. 2004). 

1.2 The Complexity of the "HB Phenomenon" 

It is virtually impossible to write a short review pa- 
per on HB stars covering "observations" , "theory" , and 
"implications for the formation of the Galaxy": each 
one of these subjects covers so much material that one 
could rather write separate review papers for each one 
of them. Moreover, a review of HB stars cannot be 
complete without looking into their progenitors and 
their progeny. The task of a reviewer of HB stars 
is accordingly a daunting one, and it is virtually im- 
possible to aim at completeness. In the present pa- 
per, while attempting to cover a broad spectrum of 
HB-rclated topics, we again hold no hope of provid- 
ing a complete review of the literature on these sub- 
jects. Recent reviews focusing on several more or 
less specific topics related to HB stars have been pro- 
vided by Cacciari (1999, 2003), Chaboyer (1999), dc 
Boer (1999), Moehler (2001, 2004), Sweigart (1997b, 
1999), Demarque (1999), Landsman (1999), Lee et 
al. (1999), Walker (2000), Green, Liebert, & Saffer 
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(2001), Cacciari & Clementini (2003), Bono (2003), Do 
Medeiros (2003), Piotto (2003), Maxted (2004a,b), Cas- 
sisi (2005), Storm (2006), Hebcr (2008), and Rood et 
al. (2008); and very instructive earlier reviews, cov- 
ering diverse astrophysical contexts, include those by 
Sweigart (1985, 1990, 1994), Philip (1994), Cox (1995), 
Dorman (1995), Smith (1995), Beers (1996), Stetson, 
VandenBerg, & Bolte (1996), Sarajedini, Chaboyer, & 
Demarque (1997), Fusi Pecci & Bellazzini (1997), Rood, 
Whitney, & D'Cruz (1997), and Rood (1998). Sim- 
ilarly, excellent sections focused on HB stars can be 
found in the reviews on the evolution of low-mass stars, 
Population II stars, globular clusters, and related top- 
ics by Renzini (1977, 1983), Iben & Renzini (1984), 
Castellani (1985, 1999), Caputo (1985, 1998), Renzini & 
Fusi Pccci (1988), Rood & Crocker (1989), Iben (1991), 
Zinn (1993a,b), D'Antona (1999), Feast (1999), Carney 
(2001), Harris (2001), and Gratton, Sneden, & Car- 
retta (2004b), among others. Other recent reviews by 
the present author on the subject of HB stars include 
Catelan (2004b, 2006, 2008a,b, 2009). 

In Figure 1 we show a visual color-magnitude dia- 
gram (CMD) for the Carina dwarf spheroidal (dSph) 
satellite of the Milky Way, with several different HB 
components indicated, including both a red clump and 
a red HB. The RR Lyme "gap" and the blue HB are 
also indicated. Similarly, Figure 2 shows CMDs for the 
Galactic globular cluster NGC 6752, in the visual (left 
panel) and in the U, U — V plane (right panel). These 
plots reveal the complexity of the blue tail phenomenon, 
with the positions of the HBA, HBB, and EHB com- 
ponents (see below) indicated, along with those of the 
Grundahl et al. (1999) and Momany et al. (2002, 2004) 
"jumps" and of possible blue HB gaps. The place where 
blue hook stars would be found, if present in the clus- 
ter (NGC 6752 actually appears to lack blue hook stars, 
according to Momany et al.), is also schematically in- 
dicated. 

In the next section, we briefly discuss each of these 
components in turn. 

2 The Different Constituents of the Horizontal 
Branch 

2.1 The red clump 

Red clumps originate from red giants that undergo the 
helium flash at the tip of the RGB while still hav- 
ing a total mass of more than ~ 1 M Q (but less than 
^2-2.5 M ; see, e.g., Fig. 2 in Girardi 1999); accord- 
ingly, they are commonly present in intermediate-age 
systems and old open clusters (e.g., Cannon 1970; Bro- 
cato, Di Carlo, & Menna 2001; Grocholski & Sarajedini 
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Fig. 1 V, B-V CMD for the Carina dSph galaxy, with 
the positions of red clump stars, red HB stars, the RR 
Lyrae "gap," and the blue HB indicated. The red clump 
is associated with the younger turnoff, at V ~ 23 mag, 
B — V ~ 0.25 mag. The other components derive from the 
old turnoff whose presence is indicated by the faint subgiant 
branch at V ~ 23.25 mag, B — V ~ 0.55 mag. Adapted from 
Monelli et al. (2003). 
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Fig. 2 V, B-V {left panel) and U, U-V (right panel) 
CMDs for NGC 6752. In the left panel, the positions of 
HBA, HBB, and EHB stars (collectively forming the "blue 
tail") are indicated, along with the positions of possible gaps 
along the distribution. In the right panel, besides the blue 
tail, the positions of the Grundahl et al. (1999) and Momany 
et al. (2002, 2004) "jumps" are also indicated, along with 
the approximate position where "blue hook" stars would be 
found if they were present in the cluster. Adapted from 
Momany et al. (2002). 



2.2 The red HB 



2002). They are physically related to, but should not 
be confused with, the so-called "red HB stars," which 
are less massive, significantly fainter and usually fall 
along the "horizontal" part of the HB, as clearly shown 
in Figure 1. 

A common mistake is to confuse "red clump" and 
"RGB bump" stars, but they are totally different types 
of stars: the latter are bona-fide red giants with a 
partially degenerate helium core that burn hydrogen 
in a thin shell surrounding the core (Thomas 1967). 
This hydrogen-burning shell moves outward in mass as 
the star evolves up the RGB, and eventually reaches a 
chemical discontinuity left behind by the previous, max- 
imum inward penetration of the hydrogen-rich convec- 
tive envelope, leading to a sudden replenishment of the 
shell with fresh nuclear fuel and a momentary reversal 
of evolutionary path on the RGB (see Salaris, Cassisi, 
& Weiss 2002 for a review and additional references). 
Recent papers discussing several different aspects of red 
clump stars include those by Girardi (1999) and Salaris, 
Percival, & Girardi (2003). 



As just stated, red HB stars are the lower-mass ana- 
logues of the red clump stars. They are commonly 
present in metal-rich (e.g., Armandroff 1988; Ortolani 
et al. 1995) as well as relatively young globular clus- 
ters (e.g., Stetson ct al. 1989; Buonanno et al. 1993). 
However, red HB stars can also be found in metal-poor 
systems with "normal" ages, where they are usually in- 
terpreted as either the progeny of red giants that lost 
very little mass on their ascent of the RGB, or as former 
blue HB stars or RR Lyrae variables that have evolved 
to the right on the color-magnitude diagram, now ap- 
proaching the end of their evolution as HB stars and 
the beginning of the asymptotic giant branch (AGB) 
phase. Another route towards producing (admittedly 
fewer) red HB stars is through the evolution of blue 
straggler stars (BSS; Renzini & Fusi Pecci 1988; Fusi 
Pecci et al. 1992). Fusi Pecci et al. suggest, in fact, 
that a red HB star "of BSS origin" should be present 
for every ~ 6 BSS present in a given globular cluster. 
(Note that, depending on the amount of mass lost by 
the BSS on the RGB phase, at least some of these might 
better classify as red clump stars.) 
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2.3 The RR Lyrae "gap" 

This is the part of the HB that crosses the Cepheid 
instability strip. The term "gap" is very inadequate, 
but is still commonly used. The reason for this is that, 
in order to properly place an RR Lyrae variable in the 
color-magnitude diagram, one needs to follow its whole 
pulsation cycle and thereby obtain reliable mean colors 
and magnitudes. Since most color-magnitude diagram 
studies lack the sufficient time coverage, these stars 
are often simply omitted from the published CMDs, 
thereby leading to an artificial "gap" at the instability 
strip level. However, there are indeed several globular 
clusters — the so-called bimodal-HB globular clusters — 
which do seem to have relatively few RR Lyrae stars 
compared to the nonvariablc HB stars to their right 
and left in the CMD (see Catelan et al. 1998 for a re- 
view and extensive references). The most famous such 
cluster is NGC 2808 (Harris 1974), which in spite of the 
recent discovery of a sizeable population of RR Lyrae 
stars, still remains firmly classified as having a bimodal 
HB (Corwin et al. 2004). 

2.4 The blue HB 

As clearly revealed by its name, blue HB stars are 
HB stars falling to the blue of the RR Lyrae insta- 
bility strip. There have been numerous subdivisions of 
the blue HB, the most common including HBA (or A- 
type HB), HBB (B-type HB), and EHB (extreme or 
extended HB) stars (see Fig. 2). HBA stars are blue 
HB stars cooler than about 12,000 K; HBB stars in- 
clude those with temperatures in the range between 
12,000 K and 20,000 K; and EHB stars include HB 
stars hotter than 20,000 K. The latter cover a remark- 
ably small range in envelope masses, generally less than 
0.02 M Q (e.g., Dorman, Rood, & O'Connell 1993)— and 
therefore also in total masses, since the He-core mass is 
essentially the same for all low-mass stars with a given 
chemical composition (e.g., Caputo & Degl'Innocenti 
1995). In a visual CMD, the blue HB component may 
contain a "horizontal" part — the canonical blue HB (see 
Fig. 1) — and an effectively "vertical" component, com- 
monly referred to as the the blue HB tail (see Fig. 2). 

Many authors have suggested that the blue HB 
proper and the blue HB tail are separated by a "gap" 
which is indeed seen in several globular clusters around 
(B — V)o ~ (e.g., Buonanno, Corsi, & Fusi Pecci 1985; 
Caloi 1999) — but perhaps not in all clusters containing 
blue tails (Catelan et al. 1998; Lee & Carney 1999a). 
On the other hand, several additional "gaps" may be 
present along the blue tail (Newell 1973; Newell & Gra- 
ham 1976; Newell & Sadler 1978; Lee & Cannon 1980; 



Ferraro et al. 1998; Piotto et al. 1999; Momany et al. 
2004); a recent, detailed description of these gaps has 
been provided by Momany et al.. Several of these gaps, 
as discussed by Catelan et al., may still require more 
sophisticated statistical analyses to establish their re- 
ality beyond reasonable doubt, particularly when well- 
observed clusters seem to lack such gaps altogether, as 
in the case of M2 (NGC 7089; Lee & Carney 1999a). 
Most impressively, Behr (2003b) has recently argued 
that the field "blue HB star" sample studied by Newell 
and Newell & Graham — and which gave rise to the very 
concept of gaps along the blue HB — is mostly com- 
prised of stars in different evolutionary phases. We 
quote Behr: 

We expected to find [based on detailed spectral analysis] 
a high fraction of HB candidates among the faint blue high- 
latitude stars listed by Newell (1973) and Newell & Graham 
(1976), especially in light of the 'gaps' in the color distribution 
of these stars... But fewer than half (11 of 27) of the Newell 
stars that we observed were clearly HB objects, with another 
11 stars classified as Population I dwarfs, and the remaining 
five stars marked as pAGB [post-AGB], subgiants, and such. 

The EHB component is the globular cluster analog of 
the field blue subdwarf (or sdB) stars (Caloi 1972). The 
majority of the (nearby) field sdB stars, according to 
Altmann, Edelmann, & de Boer (2004) and Arifyanto 
et al. (2005), appear to have disk kinematics, unlike 
HBA stars, which they find to be mostly halo stars (see 
also Altmann, Catelan, & Zoccali 2005). Accordingly, 
Altmann et al. (2004) pose the question whether field 
HBA and EHB stars truly have a similar physical origin. 

Detailed studies of the color-magnitude diagrams of 
globular clusters in several different bandpasscs indi- 
cate the presence of additional "fine structure" on the 
blue HB, most notaly the Grundahl jump (Grundahl 
et al. 1999) at a temperature around 11,500 K, which 
is seen as a sudden deviation, towards brighter magni- 
tudes, of HB stars hotter than this point, particularly 
evident when using the Stromgrcn u band (Grundahl 
et al. 1999). However, it is also present when using the 
Johnson U passband (Newell 1970; Siegel et al. 1999; 
Bedin et al. 2000; Baev, Markov, & Spassova 2001; 
Markov, Spassova, & Baev 2001; Catelan et al. 2002a; 
Momany et al. 2002, 2004), as can be clearly seen from 
Fig. 2. This "jump" has been interpreted by Grun- 
dahl et al. as due to radiative levitation and diffusion 
effects, which lead to the metal-poor HB stars hotter 
than 11,500 K possessing strongly metal- enhanced and 
helium- depleted atmospheres. 1 The earlier literature 



1 Fortunately, and as emphasized by Landsman (1999) and Behr 
(2003a), Mg appears to be virtually unaffected by radiative lev- 
itation and diffusion effects, thus allowing a "window" into the 
original mctallicity of the star — which is of extreme importance 
in the case of field stars in particular. 
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was extensively reviewed by Grundahl et al.; empiri- 
cal evidence supporting their arguments was provided 
in spectroscopic studies of blue HB stars in globular 
clusters by Bchr ct al. (1999), Moehler et al. (1999b, 
2000, 2003), Behr (2003a), and Fabbian et al. (2005). 
Diffusion calculations that showed how some of the ob- 
served spectroscopic and photometric patterns come 
about were provided by Michaud, Vauclair, & Vauclair 
(1983), Michaud, Richer, & Richard (2007, 2008), and 
Bon-Hoa, LcBlanc, & Hauschildt (2000). The "gap" 
at (B — V)o w has been tentatively associated with 
this phenomenon (Caloi 1999), but such a color appears 
much redder than would be expected for the observed 
Grundahl jump temperature of 11,500 K. In addition, 
the theoretical calculations by Grundahl et al. suggest 
that other bandpasses in the Stromgren and Johnson 
systems besides u and U should not be dramatically 
affected by this phenomenon. 

Most interestingly, recent observations have shown 
that the Grundahl jump phenomenon is also accompa- 
nied by a sharp drop in measured rotation velocities, 
blue HB stars hotter than 11,500 K presenting essen- 
tially no rotation, in contrast with cooler stars which 
may show quite significant rotation velocities, up to 
about 40kms _1 (Peterson, Tarbell, & Carney 1983; 
Peterson 1983, 1985; Peterson, Rood, & Crocker 1995; 
Behr, Cohen, & McCarthy 2000a; Behr et al. 2000b; 
Kinman et al. 2000; Rccio-Blanco ct al. 2002, 2004; 
Bchr 2003a,b; Carney et al. 2003). While Carney et 
al. note that the presence of rotation in both red and 
blue HB stars presents a problem for models in which 
rotation is a candidate second parameter, it must be 
noted that the abrupt disappearance of rotation exactly 
at the Grundahl jump temperature indicates that dif- 
fusion and levitation patterns interfere with the star's 
observed surface rotation, somehow quickly damping 
the latter for temperatures higher than about 11,500 K 
(Sills & Pinsonncault 2000; Sweigart 2002). Therefore, 
we have at present no means to check whether stars 
hotter than this temperature may have arrived on the 
ZAHB with surface rotation velocities faster than those 
observed for red HB and cooler HB stars. Asteroseis- 
mology may provide a very useful probe of the (in- 
ternal) rotation velocities of hot HB stars in the near 
future (Kawaler & Hostler 2005). A most intriguing 
piece of the puzzle is provided by the RR Lyrae stars, 
for which no evidence of rotation has been observed so 
far (see Carney et al. 2003 for extensive references). 
Why would cool blue HB stars and red HB stars show 
clear signatures of rotation, but not the intermediate- 
temperature RR Lyrae variables? More extensive spec- 
troscopic observations of RR Lyrae stars are clearly 
needed to settle this issue. 



What is the physical origin of the sudden onset of 
radiative levitation and diffusion patterns that lead to 
the Grundahl ct al. (1999) jump at 11,500 K? The an- 
swer to this question is not entirely clear at present, 
but Sweigart (2002) noted that this temperature is very 
close to the one corresponding to the disappearance of 
surface convection in HB stars, thus strongly suggesting 
a link between the disappearance of convection and the 
onset of radiative levitation and gravitational diffusion 
effects in these stars. In addition, Sweigart suggests 
that the low rotation velocities of stars hotter than the 
Grundahl jump may be due to the spin down of the 
surface layers by a weak stellar wind induced by the 
radiative levitation of Fe. As shown by Vink & Cassisi 
(2002), such winds are indeed predicted by theory. 

2.5 The extended (or extreme) HB 

Before a star arrives on the HB, it must undergo a he- 
lium flash the onset of helium burning, under partially 
degenerate conditions, which takes place at the tip of 
the RGB (Schwarzschild & Harm 1962). It is interest- 
ing to note that several early hydrodynamic studies of 
the He-flash indicated the occurrence of an explosive 
event and actual disruption of the star (e.g., Edwards 
1969; Wickett 1977; Cole & Deupree 1980). Given that 
this prediction is in clear conflict with the very exis- 
tence of HB stars, significant effort has also been de- 
voted towards the computation of hydrostatic flashes 
(e.g., Mengel & Gross 1976; Despain 1980; Mengel & 
Sweigart 1981). The hydrostatic approximation, with 
some numerical and physical refinements, is still often 
used in modern calculations (e.g., Brown et al. 2001; 
Cassisi et al. 2003; Lanz et al. 2004; Piersanti, Tor- 
nambe, & Castellani 2004; Serenelli & Weiss 2005)— 
and these have indeed been vindicated by modern hy- 
drodynamic studies (e.g., Deupree 1996; Dearborn, Lat- 
tanzio, & Eggleton 2006) which revealed that the afore- 
mentioned hydrodynamical predictions were indeed in- 
correct. It is perhaps not a conclusively settled matter 
whether any of the material that is nuclearly processed 
during the He-flash may reach the surface of the star, 
and also whether one or more mass loss episodes may 
be triggered by the primary and secondary core flashes; 
however, and to the extent that the latest He-flash cal- 
culations provide a realistic description of the process, 
no such phenomena should be expected. Indeed, Dear- 
born et al. have recently pointed out that, due to ex- 
pansion velocities that are much lower than the local 
sound speed, hydrostatic modeling should indeed cap- 
ture the essence of the helium flash process. On the 
other hand, the 3D models by Dearborn et al. intrigu- 
ingly reveal motions "of an apparently convective na- 
ture" beyond the H-burning shell, which these authors 
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claim not to be directly associated with the flash, but 
which may bring the products of H-burning to the sur- 
face of the star. 

On the other hand, the situation can become much 
more complex when the star loses so much mass on its 
ascent of the RGB that the He-flash ends up taking 
place not at the RGB tip, but rather during the helium 
white dwarf cooling curve — a so-called late hot flasher 
(Brown ct al. 2001; Cassisi ct al. 2003; Castellani, 
Castellani, & Prada Moroni 2006). In this case, exten- 
sive mixing between the envelope and regions that un- 
derwent significant hydrogen and helium burning is in- 
deed expected, and these stars, when they finally man- 
age to settle on the zero-age HB (ZAHB), may end up as 
the so-called blue hook stars (e.g., D'Cruz ct al. 1996, 
2000; Whitney et al. 1998; Rosenberg, Rccio-Blanco, 
& Garcia-Marm 2004; Busso et al. 2007; Ripepi et al. 
2007), which have temperatures in excess of 35,000 K 
(Moehler et al. 2002, 2004, 2007). This is to be com- 
pared with the case in which the star undergoes an early 
hot flash — i.e., prior to arriving on the white dwarf cool- 
ing sequence — in which case Brown et al. find that no 
mixing takes place and a "canonical" EHB star results. 
In this scenario, the reason why the blue hook stars 
appear fainter than its peers on the EHB is twofold. 
First, they tend to have smaller helium-core masses, 
due to the fact that they leave the RGB before igniting 
helium in their core, as a consequence of extreme mass 
loss on the RGB (e.g., Castellani & Castellani 1993; 
D'Cruz ct al. 1996; Brown ct al. 2001). Second, and 
most importantly in the late flasher scenario, their at- 
mospheres present large enhancements in both helium 
and carbon as a consequence of mixing, thus affecting 
the bolomctric corrections (Brown et al. 2001). 

Spectroscopic evidence largely favoring this "late 
flash mixing" scenario for the origin of blue hook stars 
has recently been presented by Lanz et al. (2004) and 
Moehler et al. (2002, 2004, 2007). An intriguing indica- 
tion from these spectroscopic studies for at least some 
of the hotter EHB stars is that some hydrogen still 
remains in their atmospheres, which is not expected 
on the basis of the theoretical models. As argued by 
Moehler ct al., some residual hydrogen may perhaps 
survive a late hot flash and later diffuse to the surface 
during the HB phase. Note also that, in analogy with 
the Grundahl et al. (1999) "jump," another "jump" 
has recently been suggested to be present in blue-tail 
globulars, namely the Momany jump (see Fig. 2), at 
a temperature around 21,000 K (Momany et al. 2002, 
2004) — which these authors conjecture to be related to 
the early helium flashers (but see Catelan 2008a for a 
critical discussion). 

Recently, a sub-population of cluster stars with a 
large, primordial enhancement in the helium abundance 



has also been suggested as a possible channel produc- 
ing the blue hook stars (Lee et al. 2005b; Yoon et al. 
2008). However, it is unclear how the abundance pat- 
terns observed in blue hook stars would be accounted 
for in this scenario (e.g., Moehler ct al. 2007). In like 
vein, Castellani et al. (2006) have recently advanced the 
intriguing suggestion that at least some of the very hot 
and undcrluminous blue hook stars might actually be 
more straightforwardly explained as photometric blends 
and/or binary stars. 

3 Variable Stars on the HB 

Until recently, it was thought that RR Lyrae stars were 
the only class of variable stars on the HB. The situation 
has changed recently, with the discovery (in a beautiful 
example of theoretical modelling preceding the observa- 
tions) of non-radial pulsation among sdB stars. These 
non-radial pulsators are now divided into two differ- 
ent groups. Therefore, there are currently three known 
classes of variable HB stars, namely: 

• RR Lyrae stars: These are radial pulsators with pe- 
riods in the range between about 0.2 d and 1.0 d, 
known to pulsate primarily in the fundamental mode 
(RR Lyrae stars of type a and b, now lumped to- 
gether as RRab or RR0 stars) and in the first over- 
tone (RR Lyrae stars of type c, nowadays also re- 
ferred to as RR1 stars) . They have long been known 
to be present in large numbers in Galactic globular 
clusters (Pickering & Bailey 1895), and were first cor- 
rectly identified as radially pulsating stars by Shapley 
(1914). The letters "a," "b," "c" were first used by 
Bailey (1902) (see his p. 132) to classify RR Lyrae 
stars as a function of their light curve shapes; Bailey 
also noted that the RRb subclass is 

... similar to (the RR)a (subclass), of which it may be 
regarded as a modification. 

Indeed, both are now known to be fundamental-mode 
pulsators; in fact, the identification of RRab's with 
fundamental pulsation and of the RRc's with first 
overtone pulsation seems to have first been clearly 
made by Schwarzschild (1940). The notation in 
which the numbers and 1 are used as opposed to 
the letters "ab" and "c," respectively, was first intro- 
duced by Alcock ct al. (2000). 

In addition, there are also double-mode pulsators 
(RRd or RR01 stars), pulsating simultaneously in the 
fundamental and first-overtone modes, as first recog- 
nized by Jerzykiewicz & Wenzel (1977) among field 
stars and by Sandage, Katem, & Sandage (1981) (see 
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their Sect. Ilia) in globular clusters (see also Gru- 
berbauer et al. 2007 for an impressively detailed fre- 
quency analysis, using MOST [Microvariability Os- 
cillations of Stars] satellite observations, of the RRd 
star AQ Leo, the prototype of this class). The nota- 
tion "RRd" for these stars appears to first have been 
used by Nemec (1984, 1985). The ratio between the 
first-overtone and fundamental period for the RRd 
stars is quite well defined; indeed, on the basis of 
the data for M68 (NGC 4590), IC 4499, and M15 
(NGC 7078) compiled in Table 8 of Kovacs & Walker 

(1999) , we obtain (Pi/P ) = 0.7454 for 38 stars, 2 
with a minimum value of Pi/Pq = 0.7433 and a max- 
imum value Pi/Po = 0.7481. For the LMC RRd, 
Alcock et al. (2000) find comparable period ratios, 
but with upper and lower values shifted downward 
by « 0.0015 (compare with their Fig. 5). An impres- 
sive summary of period ratios for a variety of stellar 
systems (the LMC, the SMC, Draco, Sculptor, M15, 
M68, IC 4499, M3, and the Galactic field) is pro- 
vided in Figure 1 of Popielski, Dziembowski, & Cas- 
sisi (2000), which basically confirms the above range 
in period ratios. Recently, Clementini et al. (2004) 
have reported on the discovery of two RRd variables 
in M3 whose period ratios, in the range 0.738-0.739, 
fall well below those for previously known RRd stars. 

It has also been suggested that RR Lyrae stars 
may pulsate in the second overtone (e.g., Demers 
& Wchlau 1977; Alcock et al. 1996; Walker & Ne- 
mec 1996; Kiss et al. 1999; Clement & Rowe 2000); 
accordingly, short-period, low-amplitude variables 
which are suspected of being second-overtone pul- 
sators are now often classified as RRe or RR2 stars 
(but see Kovacs 1998, Catelan 2004b, Bono et al. 
1997a; 1997b, for arguments suggesting that at least 
some of them may rather represent the short-period 
end of the RRc distribution). In fact, Alcock et al. 

(2000) suggest that double-mode variables pulsating 
simultaneously in the first and second overtones may 
also exist, and tentatively assign them an RR12 sub- 
class. 

Non-radial modes have also been suggested to be 
present in a fraction of the RR Lyrae stars (e.g., 
Kovacs 1995; Kolcnberg 2002; Gruberbauer et al. 
2007; and references therein), primarily as a means 
to explain the so-called Blazhko (1907) effect (e.g., 
Dziembowski & Mizerski 2004). The Blazhko effect 
consists in a periodic modulation, on a much longer 
timescale than the primary period, of the light curve 



2 It is basically this ratio that allows one to compute the "funda- 
mentalized" period of an RRc star by using the relation log Po = 
log Pi +0.128. 
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shape. 3 The modulation (or Blazhko) period falls 
in the range between 5.309 d (as found by Jurcsik 
et al. 2006 for the field star SS Cnc) and 530 d (as 
found by Nagy 1998 for the field star RS Boo). A 
variety of other multiperiodic phenomena may also 
take place in RR Lyrae stars (sec, e.g., Alcock et al. 
2000, 2004; Mizerski 2004; Gruberbauer et al. 2007). 
We note, in passing, that the Fourier decomposition 
(D m ) method (Jurcsik & Kovacs 1996) has recently 
been criticized by Cacciari, Corwin, & Carney (2005) 
as a diagnostic of the Blazhko effect in RR Lyrae 
stars. 

• V361 Hya variables: Also known as EC 14026 or 
sdBV stars, these are non-radial, short-period, low- 
order p-modc pulsators, whose periods fall in the 
range 80-400 s, and whose amplitudes bracket the 
interval between 4 and 25 mmag. Their tempera- 
tures are found in the range between 29,000 K and 
36,000 K, and their gravities in the range 5.2 < 
log g < 6.1 (Kilkenny 2002; Fontaine et al. 2004, 
2006a, b). Note that this class of variable star had 
been predicted (Charpinet et al. 1996) before it was 
first observed (Kilkenny et al. 1997). 

• V1093 Her variables: Also known as PG1716+426 
or "Betsy" variables (after Elizabeth Green, the dis- 
coverer of the group), these are non-radial, relatively 
long-period, high-order g-mode pulsators (e.g., Green 
et al. 2003; Fontaine et al. 2003; Reed et al. 2004). 
Their periods fall in the range 2000-9000 s, and their 
amplitudes are very small, generally being smaller 
than 0.5 mmag. Their temperatures bracket the in- 
terval 25,000-30,000 K, and their gravities are in the 
range 5.1 < log g < 5.8 (Fontaine et al. 2004, 2006a,b; 
Reed et al. 2004). 

• Mixed-mode variables: Variable stars showing both 
V361 Hya and V1093 Her characteristics have re- 
cently been discovered as well (Baran et al. 2005; 
Schuh et al. 2006). 

For a recent discussion of the pulsation mechanism 
(the "Fe-bump opacity mechanism" ) driving the pulsa- 
tions in V361 Hya and V1093 Her stars, see Jeffery & 
Saio (2006). 

It is important to emphasize that the presence of 
variable stars on the HB phase provides us with a 
unique opportunity to utilize stellar pulsation obser- 
vations and theory to improve our understanding of 



3 In their review of 100 years of observations of the star RR Lyrae, 
Szeidl & Kollath (2000) point out that the phenomenon might 
more appropriately be called the Blazhko- Shapley effect, since 
Shaplcy (1916) was actually the first to demonstrate that the 
oscillations in the maxima of an RR Lyrae star can be described 
as a periodic variation in the shape of the light curve and in the 
height of the maxima. 
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HB stars in general. As well known, the pulsation 
properties of stars are fundamentally related to their 
mean densities through the so-called period-mean den- 
sity relation — an expression originally due to Rittcr 
(1879). This extremely important equation specifies 
that the period is inversely proportional to the square 
root of the mean density of the star, which in turn can 
be expressed in terms of the star's global physical pa- 
rameters, such as mass M, luminosity L, and effective 
temperature T e g (e.g., van Albada & Baker 1971; Bono 
et al. 1997b; Caputo, Marconi, & Santolamazza 1998; 
Di Criscienzo, Marconi, & Caputo 2004). As a mat- 
ter of fact, Ritter's relation is applicable, to first order, 
over an impressively wide range of stellar parameters 
(e.g., Cox 1974). In addition, the observed periods — 
and especially so in the case of non-radial pulsators — 
depend on the star's detailed structural profile. This 
is a natural consequence of the fact that the measured 
periods are directly related to the speed of travel of the 
sound waves across the stellar interior. In this sense, 
the importance of the non-radial variables on the HB to 
constrain the internal structures of HB stars, including 
the amount of internal rotation, has recently been em- 
phasized by Kawaler (2006), Kawaler & Hostler (2005), 
and Fontaine et al. (2006a), and nicely demonstrated 
in practice by Charpinet et al. (2005, 2006), Randall et 
al. (2007), and Van Grootel et al. (2008). It is inter- 
esting to note that the high rotation velocity recently 
derived by Dixon, Landsman, & Brown (2004) for the 
post-AGB star ZNG 1 in M5 might be accounted for, as 
noted by those authors, if RGB and HB stars are able to 
maintain rapid rotating cores, as may also be required 
(Sills & Pinsonneault 2000) to explain the presence of 
fast rotators on the HB (see Sect. 2.4 for extensive ref- 
erences) . 

It is worth noting that several authors have called 
attention to the possible presence of non-variable 
HB stars inside the RR Lyrae instability strip (e.g., 
Sandagc & Katcm 1968; Wehlau 1990; Xiong, Cheng, 
& Deng 1998; Stetson, Catelan, & Smith 2005; and ref- 
erences therein) . The recent discovery of Cepheid vari- 
ables presenting pulsation amplitudes at the mmag level 
(Buchler et al. 2005) suggests that at least some of these 
"non-variable" stars may indeed be varying, though 
with much smaller amplitudes than would have been 
possible to detect with more traditional techniques — as 
nicely demonstrated by Clement & Rowe (2001) in the 
case of NGC 5897. According to Buchler et al., "ultra- 
low amplitude" RR Lyrae stars are predicted to exist 
close to both the blue and red edges of the instability 
strip. Note, on the other hand, that the General Cat- 
alog of Variable Stars (Kholopov et al. 1998) defines 
RR Lyrae stars as having amplitudes in the range 0.2 



to 2 mag in V, 4 thus showing that RR Lyrae with V 
amplitudes smaller than 0.2 mag, though known to ex- 
ist even among the normally higher-amplitude RRab's 
(e.g., Sandage 1981; Grenon & Waelkens 1986; Clement 
& Rowe 2001; Soszyhski et al. 2002, 2003; Cacciari et 
al. 2005), are indeed quite rare. Note that variable stars 
falling on the red HB — i.e., outside the formal bound- 
aries of the RR Lyrae instability strip — have also been 
suggested to exist (Xiong, Cheng, & Deng 1998 and 
references therein). 

Still in regard to pulsation amplitudes, an interest- 
ing recent result worth mentioning are the very large 
amplitudes — up to ~ 5 mag — found for RR Lyrae stars 
in the far UV (Downes et al. 2004; Dieball et al. 2005, 
2007; Welsh et al. 2005; Wheatley et al. 2005; see also 
Fig. 1 in Bonnell et al. 1982 for an example of an ear- 
lier RR Lyrae light curve in the far UV, already indica- 
tive of large amplitudes), and up to <~ 2.6 mag in the 
near UV (Browne 2005). Conversely, in the near-IR, 
RR Lyrae amplitudes become quite small (e.g., Jones, 
Carney & Latham 1998; Liu & Janes 1990a). Thus, 
RR Lyrae stars usually reveal themselves more eas- 
ily when time-series surveys are conducted using bluer 
bandpasses, whereas average magnitudes can be more 
easily computed using near-IR observations (see also 
§6.3.1 below). 

Before closing this section, we would like to com- 
ment on the possibility that some stars in the rapid 
phase of evolution between the RGB tip and the ZAHB 
do in fact become variables. Evolutionary paths for 
these stars are shown, for instance, in Fig. 4 in Mengel 
& Sweigart (1981) and Fig. 4 in Brown et al. (2001), 
where it can be seen that many of these stars do cross 
the instability strip before settling on the ZAHB. While 
distinguishing them from more ordinary RGB, AGB, 
and RR Lyrae stars is obviously far from straightfor- 
ward, both due to the very small expected number of 
stars in this phase and because of their overlapping with 
them on the CMD, their extremely fast evolutionary 
timescale (~ 10 6 yr) hints at the possibility of identify- 
ing them by appropriate monitoring in wide-field sur- 
veys. This is because much more extreme period change 
rates might obtain for them than for most other vari- 
ables in a similar position on the CMD. To the best of 
our knowledge, no other technique has been proposed 
in the literature for the detection of these elusive post- 
RGB tip/pre-ZAHB stars. The reader is referred to 
the recent paper by Silva Aguirre et al. (2008) for an 
application of this method to the case of the Galactic 
globular cluster M3. 



Sec http://www.sai.rasu.su/groups/cluster/gcvs/gcvs/iii/vartype.txt. 
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4 Detecting HB Stars in Distant Systems: 
Methods 

While HB stars are relatively bright, their direct de- 
tection in distant systems (i.e., at the distance of M31 
and beyond) is a challenge. This is particularly so in 
the case of blue HB stars, which can be several magni- 
tudes fainter than the "horizontal" level in V (Fig. 2), 
due to the increase in the U-band bolomctric correc- 
tion at high temperatures. While the blue HB stars be- 
come brighter towards the far ultraviolet, which could 
make them candidates for detection from space, the 
required exposure times, at the distance of M31, are 
prohibitively large. Therefore, while extremely long- 
exposure HST observations have demonstrated that it 
is possible to reliably detect blue HB stars down to 
at least the main sequence turnoff level in M31 (e.g., 
Brown et al. 2003, 2004b, 2006), alternative approaches 
are needed at present to obtain a more complete assess- 
ment of the HB morphologies in extragalactic systems. 

4.1 Far-UV Observations 

The far-UV output from globular clusters is strongly 
dependent on the temperature of the sources. Since 
hot HB stars are important UV emitters (e.g., Dor- 
man, O'Connell, & Rood 1995; see also Catelan 2008b 
for a recent review), this is expected to translate into 
a far UV (integrated) color-HB type correlation, espe- 
cially when the contribution of bright, individual far- 
UV sources (such as post-AGB stars) , which are present 
in clusters in non-statistically significant numbers, is re- 
moved from estimates of the integrated far UV colors. 

Landsman et al. (2001) and Landsman & Catelan 
(2009, in preparation) tested this prediction, using in- 
tegrated far UV fluxes for Galactic globular clusters as 
summarized by Dorman et al. (1995), and including re- 
visions to the far UV photometry based on images taken 
with the Ultraviolet Imaging Telescope (UIT). The re- 
sult is shown in Figure 3, which depicts the trend of 
variation in the (15 — V)o color (where "15" stands for 
a bandpass centered at 1500 A) as a function of the 
HB morphology parameter (B2-TZ)/(B+V+1Z) (Buo- 
nanno 1993; Buonanno et al. 1997), where B2 is the 
number of blue HB stars bluer than (B — V)o — —0.02, 
and £>, V, 1Z arc the numbers of blue, variable (RR 
Lyrae), and red HB stars, respectively. The latter 
number counts are also widely used in constructing the 
so-called Lee-Zinn parameter, C — (B—1Z)/(B+V+1l) 
(Zinn 1986; Lee 1990; Lee et al. 1990). The reason why 
we use the Buonanno as opposed to the Lee-Zinn pa- 
rameter is that it is more successful at breaking the 
degeneracy that characterizes C in the case of clus- 
ters with completely blue HBs. More specifically, C is 
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(B2-R)/(B+V+R) 

Fig. 3 Variation in the integrated (15— V)o color of Galac- 
tic globular clusters as a function of the HB morphol- 
ogy parameter (B2-TZ)/(B+V+Tl) (Buonanno 1993). The 
shaded area (labeled "bimodal HB") indicates a region oc- 
cupied by clusters which have sizeable blue HB components, 
but which at the same time have sizeable red HB compo- 
nents as well. The lines indicate theoretical predictions 
based on solar-scaled (continuous lines) and a-enhanced 
(dotted gray lines) abundances. From Landsman et al. 
(2001) and Landsman & Catelan (2009, in preparation). 

not able to distinguish between a long blue tail cluster 
such as M13 or NGC 6752 and a stubby blue HB clus- 
ter such as NGC 288 (all having C ~ +1.0), whereas 
(B2-K)/(B+V+K) (as well as several other HB mor- 
phology parameters, including some of those defined 
by Fusi Pecci et al. 1993 and Catelan et al. 2001a) arc 
better suited for the task. Indeed, Figure 3 shows that, 
apart from the clusters with bimodal HBs (NGC 6388, 
NGC 6441, M75), there is a good correlation between 
the Buonanno parameter and (15 — V)o, characterized 
by a high correlation coefficient r = 0.93. Using C, 
one finds an r = 0.82 instead, but this is primarily due 
to the large difference in far UV color between 47 Tu- 
canae (NGC 104) and the group of blue HB clusters: 
removing 47 Tuc, the correlation coefficient, when using 
C, drops to r = 0.57, whereas the one obtained using 
the Buonanno parameters still remains fairly high at 
r = 0.84. These correlations are in good agreement 
with the theoretical predictions by Landsman et al. for 
an a-enhanced composition, as can clearly be seen from 
the dotted line in Figure 3. In summary, one should 
be able to estimate, at least as a first approximation, 
the HB type of an old extragalactic globular cluster by 
measuring its integrated far UV light. 

Indeed, that the far UV light is a powerful diagnos- 
tic of the presence of hot HB stars has been demon- 
strated by the observation of a peculiarly high far UV 
flux from the moderately metal-rich Galactic globular 
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clusters NGC 6388 and NGC 6441 (Rich, Minniti, & 
Liebert 1993), which led the authors to suggest the ex- 
istence of long blue tails in these clusters several years 
before these were first directly observed with HST (Pi- 
otto et al. 1997; Rich et al. 1997). Recently, Peterson et 
al. (2003) have performed a detailed analysis of the mid- 
UV flux from the extremely massive M31 globular clus- 
ter Gl, finding that it likely presents both regular and 
extreme blue HB stars. Its recent detection in far-UV 
observations with GALEX (Rey et al. 2005, 2007) lends 
support to their conclusion. The HST color-magnitude 
diagram for the cluster (Rich et al. 2005) cannot reveal 
the presence of hot HB stars since it just reaches the HB 
level, but it is not clearly inconsistent with the presence 
of a few blue HB stars at the "horizontal" level. 5 

To close, it is important to note that the H/3 and 
US integrated indices can be sensitive indicators of the 
presence and temperature of HB stars in unresolved 
stellar systems (e.g., Lee, Yoon, & Lee 2002; Maras- 
ton et al. 2003; Schiavon et al. 2004; Percival et al. 
2009; and references therein). The dependence of these 
and several other photometric and spectroscopic indi- 
cators on HB morphology has recently also been dis- 
cussed by Lee, Lee, & Gibson (2002), Proctor, Forbes, 
& Beasley (2004), Salaris & Cassisi (2007), and Koleva 
et al. (2008), among others. 

4.2 Type II Cepheids: Bright Indicators of the 
Presence of Faint Blue HB Stars 

Another technique to infer the presence of a blue HB 
component in extragalactic systems without the need 
to obtain exceedingly deep photometry could use fairly 
bright type II Cepheids (including BL Herculis, W Vir- 
ginis, and RV Tauri stars) as an indicator. Indeed, 
type II Cepheids are widely believed to be the immedi- 
ate progeny of HB stars with little envelope mass (but 
with still enough to reach the AGB stage) — i.e., blue 
HB stars. This is supported both by the observational 
record, which shows that type II Cepheids are present 
only when a sizeable blue HB component is also present 
(Wallcrstcin 1970; Smith & Wehlau 1985), seemingly ir- 
respective of the metallicity (Pritzl et al. 2002b, 2003) ; 6 
and by theoretical models (Schwarzschild & Harm 1970; 
Gingold 1976; 1985; see also Di Criscicnzo et al. 2007 



5 Whilc Rich ct al. (2005) provide C values for their observed 
M31 globular clusters, no completeness corrections were applied 
in deriving them, which may lead to important underestimates 
of the number of blue HB stars — and therefore of C — for those 
among them possessing well-populated blue HB tails. 

6 There is a single known exception to this rule, provided by 
the type II Cepheid in Palomar 3 (Borissova, Ivanov, & Cate- 
lan 2000). 



for a recent discussion and additional references). Ac- 
cordingly, detection of type II Cepheids should imme- 
diately imply the presence of a sizeable blue HB com- 
ponent. In this sense, it is worth noting that Dolphin 
et al. (2004) and Williams (2005) found possible type II 
Cepheids in variability studies of small M31 halo fields; 
if this classification is confirmed (the authors note that 
the stars could also be Anomalous Cepheids, for in- 
stance), this would be consistent with the deep M31 
CMDs obtained by Brown ct al. (2003, 2006), which 
clearly reveal the presence of a sizeable (though seem- 
ingly not very extended) blue HB component. Like- 
wise, Pritzl et al. (2005a) suggest that the M31 dSph 
satellites And I and And III may also contain a small 
number of type II Cepheids, which should also imply 
the presence of blue HB components. Their HST CMDs 
for these galaxies do indeed suggest the presence of blue 
HB stars, particularly evident in the case of And I. Last 
but not least, we note that Fliri et al. (2006) and Vi- 
lardel, Jordi, & Ribas (2007) have recently reported the 
detection of type II Cepheids in M31, including many 
in the direction of the M31 bulge — which, if confirmed, 
implies the presence of a sizeable blue HB component 
in that direction as well. 



5 Mass Loss on the RGB 

Knowledge of how red giants lose mass is one of the 
indispensable conditions for understanding the proper- 
ties of HB stars, including their color distributions, the 
production of RR Lyrae and sdB stars, the variation in 
HB morphology with metallicity, and the (in)famous 
second-parameter phenomenon. Several previous re- 
viewers have emphasized the need for an advancement 
in our knowledge of RGB mass loss for a consensus 
on a variety of problems involving HB morphology to 
be achieved (e.g., Rood 1973, 1998; Rood & Crocker 
1989; Rood et al. 1998). Unfortunately, even though 
more than three decades have passed since a very con- 
venient mass loss formula was advanced by Reimers 
(1975a, b), his expression is still widely used as a "law," 
even though it is now known that the Reimers mass 
loss formula does not properly account for the mass 
loss rates that have become available in the more re- 
cent literature. In fact, several other formulae have 
been proposed that do provide a better description of 
these empirical mass loss rates, some of which are sum- 
marized in Table 1. Additional formulations and ref- 
erences are provided in Origlia et al. (2007), Schroder 

6 Cuntz (2007), and McDonald & van Loon (2007). 
Boyer et al. (2008) provide a critical discussion of re- 
cent results based on Spitzer Space Telescope infrared 
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Formula name Formula for dM ^ t GB (in M Q yr 1 ) 
Reimers 5.5 x KT 13 

Modified Reimers 8.5 x KT 10 (j^j 



, +1.4 



Mullan 2.4 x KT 11 (^72) 
Goldberg 1.2 x 10" 15 i?+ 3 ' 2 

Judge-Stencel 6.3 x KT 8 g~ 1Ji 

VandenBerg 3.4 x KT 12 L +11 g-°- 9 



"Reimers formula as given in Kudritzki & Reimers (1978). All other expressions derived by Catelan (2000see his Appendix A for 
details) on the basis of the data provided by Judge & Stencel (1991). The gravity g is in cgs units, and luminosity L and radius R in 
solar units. 
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Fig. 4 Variation in integrated mass loss along the RGB 
with metallicity, for an age of 12 Gyr and several different 
prescriptions for mass loss (see Table 1). For all equations, 
mass loss values were normalized to a value of AMrgb = 
0.10 M at [Fe/H] = -1.5. 

observations of Galactic globular clusters, and caution 
that some of the results by Origlia et al. may actually 
be due to blends in their infrared data. 

Figure 4 shows the variation in the integrated mass 
loss along the RGB, on the basis of different prescrip- 
tion for the mass loss rate, for an age of 12 Gyr and 
several different metallicities (see Sect. 10 below for 
more information on the RGB models used in these 
calculations, and Catelan 2009 for an extension to the 
Origlia et al. 2007 and Schroder & Cuntz 2007 cases). 



Clearly, the metallicity dependence differs substantially 
among the different expressions. This can have impor- 
tant consequences for the prediction of the temperature 
distribution of hot HB stars in metal-rich systems, with 
implications not only for our understanding of the ori- 
gin and nature of hot HB stars in the Galactic bulge 
(e.g., Busso et al. 2005; Busso & Moehler 2008), but also 
for our modelling and interpretation of the UV-upturn 
phenomenon in elliptical galaxies and bulges of spirals 
(e.g., Horch, Demarquc, & Pinsonneault 1992; Bressan, 
Chiosi, & Fagotto 1994; Dorman et al. 1995; Yi, De- 
marque, & Kim 1997; Yi et al. 1999), as well as for our 
understanding of the relatively red HB morphologies of 
the most metal-poor globular clusters in our galaxy (see 
Fig. 7 below). In Sect. 7, we shall also address the pos- 
sible impact of several different mass loss formulae upon 
analyses of specific pairs of second-parameter globular 
clusters. 

On the other hand, the reader should bear in mind 
that it remains very unclear at present whether any 
analytical formula can be reliably used to compute the 
mass loss in red giant stars. For instance, Origlia et al. 
(2002) find, using ISOCAM near-IR data for six globu- 
lar clusters with a range in metallicities, that: i) None of 
the formulae given in Table 1 reproduces their derived 
mass loss rates, the latter being in general larger than 
the former by at least one order of magnitude; ii) There 
is no clear dependence between mass loss rate and any 
of the basic stellar parameters L, g, or R; iii) Mass 
loss takes place near the RGB tip only (more specifi- 
cally, within the very final 10 6 years of evolution on the 
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RGB — compare with the timescale shown in Fig. 16 be- 
low), and is not constant but rather episodic; iv) The 
mass loss episodes must last longer than a few days, 
but less than 10 6 yr; v) There is no clear correlation 
between mass loss and metallicity. We will discuss the 
implications of their results for the evolutionary inter- 
pretation of the second parameter problem in Sect. 7; 
see also Sect. 6.2 below for a discussion of the impact of 
different mass loss recipes upon theoretical isochrones 
in the C - [Fc/H] plane. 

Assuming their results are able to withstand the test 
of time (see Boyer et al. 2008 for some recently raised 
caveats), how does one reconcile the Origlia et al. (2002) 
study with the empirical results upon which the formu- 
lae in Table 1 are based? The answer to this question 
is not clear at present, but it is possible that the latter 
provide a better description of mass loss on the AGB, 
the former being instead more suitable for first-ascent 
giants. The connection between mass loss and stellar 
variability, which has frequently been addressed in the 
literature (e.g., Willson & Bowen 1984; Bowcn 1988; 
Willson 1988, 2000; Ramdani & Jorisscn 2001; Lebzel- 
ter & Wood 2005), may provide the key to the riddle. 
Lebzelter et al. (2005) concluded that, at a similar lu- 
minosity, red giants with higher pulsation amplitudes 
present higher mass loss rates. At the same time, while 
variability in AGB stars (which become more easily en- 
shrouded in dust) has long been established, that in 
first-ascent red giants is a fairly recent result (Ita et al. 
2002). Importantly, it appears rather clear now that 
AGB and RGB stars have different pulsation proper- 
ties (e.g., Kiss & Bedding 2004; Soszynski et al. 2004), 
first-ascent variables having smaller pulsation ampli- 
tudes and falling on a different period-luminosity rela- 
tion than AGB variables. On the other hand, Soszynski 
et al. suggest, based on the coincidence of the RGB 
pulsators belonging to the LMC, the SMC, and the 
Galactic bulge in a Petersen (1973) (or period ratio 
vs. period) diagram, that RGB pulsators with different 
metallicities do not have significantly different pulsa- 
tion properties. If confirmed, these results may provide 
very important constraints on the extent to which mass 
loss on the RGB may vary with metallicity. 

6 The Oosterhoff Dichotomy: Constraints on 
the Galaxy's Formation History 

6.1 The Oosterhoff Dichotomy: Systematics 

It has recently been argued that the Oosterhoff (1939, 
1944) dichotomy may very well hold the key to the 
formation history of the Galactic halo. The argument 



goes as follows: the Galactic halo shows a sharp divi- 
sion between Oosterhoff type I (Ool, average periods 
of the ab-type RR Lyrae variables (P a b) ~ 0.55 d) 
and Oosterhoff type II (OoII, with (P ab ) w 0.65 d) 
globular clusters, with very few clusters in the range 
0.58 < (Pab(d)) < 0.62 (the "Oosterhoff gap" ). 7 On the 
other hand, the dwarf spheroidal (dSph) satellite galax- 
ies of the Milky Way, as well as their respective globular 
clusters, fall preferentially on the "Oosterhoff gap" re- 
gion. One of the main scenarios for the formation of the 
Galactic halo envisages the build-up of the halo from 
the accretion of smaller "protogalactic fragments" not 
unlike the present-day Milky Way dSph satellite galax- 
ies (e.g., Searle & Zinn 1978; Zinn 1993b). However, if 
this were the case, the present-day halo should not dis- 
play the Oosterhoff dichotomy, since the dSph galaxies 
and their globular clusters are predominantly interme- 
diate between the two Oosterhoff classes. Therefore, 
the Galactic halo cannot have been assembled by the 
accretion of dwarf galaxies resembling the present-day 
Milky Way satellites — including, in fact, the LMC dlrr 
(Catelan 2004b). 

One criticism that might perhaps be drawn against 
this argument is related to the fact that not too many 
globular clusters satisfied the fairly strict selection cri- 
teria established by Catelan (2004b), which restricted 
his sample to globular clusters containing at least 10 
known RRab variables with measured periods. With 
this selection criterion, Catelan found that 19 globu- 
lar clusters were Ool, 9 were OoII, and two (or 6.7%) 
were Oosterhoff-intermediate. Two additional globular 
clusters, NGC 6388 and NGC 6441, were tentatively 
assigned to a new class, OoIII (Pritzl et al. 2000, 2001, 
2002b, 2003). As one can easily see, only a relatively 
small fraction of the Galactic globular clusters, of order 
20%, were included. 

In order to improve the statistics, we add to the sam- 
ple originally studied by Catelan (2004b) all Galactic 
globular clusters with at least 5 known RRab variables 



7 The question whether the Galactic halo field also presents the 
Oosterhoff dichotomy, as found by Suntzcff, Kinman, & Kraft 
(1991), or not, as suggested by the QUEST (Vivas & Zinn 2003) 
and ROTSE (Kinemuchi 2004) surveys (sec also Sandagc 2006), 
remains an open issue at present (see also Catelan 2004b). We 
note, however, that inspection of Figure 9 in Vivas et al. (2004) 
reveals that many of the candidate Oosterhoff-intermediate stars 
in the QUEST survey (e.g., their stars 1, 9, 27, 95) have ex- 
tremely uncertain amplitudes, thus rendering their position in 
the Bailey diagram — and thus their Oosterhoff status — similarly 
uncertain (Catelan 2006). The recently released results of the 
LONEOS-I survey also strongly support the reality of the Oost- 
erhoff dichotomy among field halo stars: in this sense, Figures 19 
and 20 in Miceli et al. (2008) constitute especially striking evi- 
dence. 
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with measured periods. The result is shown in Fig- 
ure 5 (left panel), where clusters with 5 to 10 RRab's 
are shown with smaller symbols than those with 10 
or more fundamental-mode variables. While the main 
source is the extensive compilation and catalog of vari- 
able stars in Galactic globular clusters by Clement et al. 
(2001), 8 this plot includes several updates compared to 
similar previous plots, such as the new measurements 
by Szekely et al. (2006) for NGC 362, by Contreras et 
al. (2005) for NGC 6266 (M62), and by Corwin et al. 
(2006) for NGC 6388 and NGC 6441. 

In addition, we separate the clusters into "bulge/ 
disk," "young halo" and "old halo" subsystems, fol- 
lowing the classification scheme by Mackey & van den 
Bergh (2005) (based on the positions of globular clus- 
ters in the HB morphology-metallicity plane). Except 
for the peculiar positions of the clusters labeled as OoIII 
in this plot, it is very clear that both the "young" and 
the "old" halo components present the Oosterhoff di- 
chotomy, perhaps the only apparent systematic differ- 
ence between the two being a systematic shift of the 
"young" clusters towards metallicities that are lower 
than for the "old" clusters, by w 0.25 dex. Impor- 
tantly, the sample size is now significantly larger, thus 
clearly improving the statistics: we find that, out of 
a total of 41 Galactic globulars with periods measured 
for at least 5 RRab stars, only 4 (or less than 10%) arc 
Oostcrhoff-intcrmcdiate. 

In Figure 5 (right panel), dwarf satellite galaxies of 
the Milky Way and their respective globular cluster 
systems are added to the previous plot. Shown are 
not only the LMC, SMC, Sagittarius dSph and Fornax 
dSph systems (along with their globular clusters), but 
also the recently discovered Canis Major dSph and its 
suggested globular cluster system (Martin et al. 2004) — 
the galaxy itself (if it indeed exists) appears to be ex- 
tremely poor in RR Lyrae variables (Kinman, Saha, 
& Pier 2004) — as well as ui Centauri. The latter, ac- 
cording to many authors, is also the remnant of a dwarf 
galaxy that was accreted by the main body of the Milky 
Way long ago (e.g., Norris et al. 1996; Dincscu 2002; 
Altmann et al. 2005; Meza et al. 2005; Da Costa & 
Coleman 2008; and references therein). This plot dra- 
matically illustrates the contrast between the Ooster- 
hoff behavior of Galactic globular clusters and the sys- 
tems associated with the dwarf satellite galaxies of the 
Milky Way. Note that the "external" systems shown in 
the figure fall more naturally along the trend defined by 
the "young halo" globular clusters — which, as already 
pointed out, are shifted with respect to the "old halo" 
globulars by about 0.25 dex, towards lower [Fe/H]. 



http: //www. astro .utoronto . ca/~cclement/read .html 



The globular cluster data used in Figure 5 are given 
in Table 2 (Galactic globulars) and Table 3 (globular 
clusters associated with the dwarf galaxy satellites of 
the Milky Way). The information contained therein 
comes from the following sources: 

• For the Fornax dSph globular clusters 3, 4, and 5, we 
adopt the (P a b) values from Greco et al. (2005; see 
also Greco et al. 2007 for a recent, detailed discussion 
of the case of Fornax 4). For clusters 1 and 2, we 
adopt the values from the entry "RRab with good P" 
in Table 4 of Mackey & Gilmore (2003). Metallicity 
values are the ones provided by Mackey & Gilmore 
(2004b). 

• For globular clusters that have been associated with 
the Sagittarius dSph and which have at least 5 
RR Lyrae variables, we have adopted (P a b) val- 
ues from Cacciari, Bcllazzini, & Colucci (2002) 
(M54 = NGC 6715), Salinas et al. (2005) (Arp 2, 
NGC 5634), and Stetson et al. (2005) (NGC 4147). 
While the association of M54 and Arp 2 to the Sagit- 
tarius dSph dates back to early studies of the sys- 
tem (e.g., Ibata et al. 1995; Da Costa & Armandroff 
1995), that of the other two quoted globular clusters 
has only recently been advanced (NGC 4147: Bellazz- 
ini, Ferraro, & Ibata 2003a; Bellazzini et al. 2003b; 
NGC 5634: Bcllazzini, Ferraro, & Ibata 2002). 9 
NGC 5634 is an interesting case, since the prelim- 
inary results of Salinas et al. (2005) indicate that, in 



9 Note that the possibility has also been raised that the massive 
outer-halo globular cluster NGC 2419 has once been associated 
with the Sagittarius dSph (Newberg et al. 2003), or perhaps been 
the stripped nucleus of a dwarf galaxy (van den Bergh & Mackey 
2004) — but we consider it a bona fide Galactic globular cluster in 
the present analysis (see also Ripepi et al. 2007). There are a few 
additional clusters whose association with Sagittarius has been 
proposed (e.g., Dinescu et al. 2001; Palma, Majewski, & Johnston 
2002; Majewski et al. 2003; Carraro, Zinn, & Moni Bidin 2007). 
However, it is well known that the Fornax dSph, with its five 
globular clusters, has an anomalously high globular cluster spe- 
cific frequency (e.g., van den Bergh 1998 and references therein); 
accordingly, Sagittarius, with a total luminosity that seems com- 
parable to that of Fornax (van den Bergh 2000; Majewski et al. 
2003), would have an even more anomalous specific frequency of 
globular clusters if all these candidates (in addition to the five 
fairly clear-cut cases of M54, Arp 2, Ter 7, Ter 8, and Pal 12) 
were indeed associated to it. Therefore, caution recommends not 
to attribute final membership status to all the new candidates 
before their association to the galaxy has been conclusively es- 
tablished. Likewise, we do not include NGC 6388 and NGC 6441 
as extragalactic globular clusters, in spite of their suggested as- 
sociation with dwarf galaxies that were long ago captured by the 
Milky Way (Ree et al. 2002). Finally, we note that Lee, Gim, & 
Casetti-Dinescu (2007) recently raised the intriguing possibility 
that pretty much all globular clusters with extended HBs are of 
extragalactic origin (being the former nuclei of dwarf galaxies), 
unlike most other globulars in the halo; the latter would be more 
consistent with a dissipational collapse scenario. 
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Fig. 5 (Left panel) The Oosterhoff dichotomy among Galactic globular clusters. Clusters belonging to the bulge or disk 
are shown as squares; those belonging to the "young halo" in the Mackey & van den Bergh (2005) classification scheme 
are shown as filled gray circles; and clusters belonging to their "old halo" are shown as filled black circles. Clusters with 
at least 10 RRab stars with measured periods are shown as large symbols, while those with 5 to 9 RRab's are shown with 
smaller symbols. (Right panel) Same as in the previous plot, but now including dwarf satellite galaxies of the Milky Way 
and their associated globular clusters. The dSph galaxies and the Magellanic Clouds are shown as open inverted triangles; 
LMC globulars are shown as open triangles; globular clusters associated with the Fornax dSph are shown as open losanges; 
those associated with the Sagittarius dSph are shown as open stars; possible Canis Major dSph globulars are shown as 
black circles with gray contours; and u> Cen is shown as an open hexagon. It is obvious from this plot that the dwarf 
galaxies orbiting the Milky Way and their associated globular clusters preferentially occupy the Oosterhoff gap region, in 
stark contrast with the Galactic globular clusters. 



spite of a (P a b) value indicative of OoII status, its 
Bailey (or period-amplitude) diagram 10 suggests in- 
stead that the cluster is Oosterhoff-intermediate (see 
Catelan 2004b for a discussion of the importance of 
the Bailey diagram in defining Oosterhoff status). 
Metallicity values are from the Harris (1996) cata- 
log 11 (Feb. 2003 update). 
• For globular clusters which have been associated with 
the CMa dSph (Frinchaboy et al. 2004; Martin ct 
al. 2004), the data come from Corwin et al. (2004) 



"The widespread usage of the term "Bailey diagram" to refer 
to the period-amplitude diagram appears to be relatively recent: 
for instance, it is not used either in Sandage's (1958) Vatican 
Conference paper or in Smith's (1995) monograph. To be sure, 
Bailey (1913, 1919) did produce plots showing, among many other 
quantities, "range" (i.e., amplitude) vs. period (his Figs. 4 and 7, 
respectively), but no great emphasis appears to have been placed 
by him on this specific diagram, other than noting that "there 
appears to be... a fairly well marked relation between the period 
and the maximum magnitude, the mimimum magnitude, and the 
range of variation" (Bailey 1913, his p. 87). 

^http : //www .physics .mcmaster . ca/Globular .html 



(NGC 2808), Walker (1998) (NGC 1851, with the ad- 
dition of five new confirmed RRab's from Sumerel et 
al. 2004), and Zorotovic et al. (2009, in preparation) 
(NGC 5286). Metallicity values are from the Harris 
(1996) catalog (Feb. 2003 update). 
• As to the LMC globular clusters, we retrieved OGLE 
data (Soszyhski et al. 2003) from their online cata- 
log, 12 and derived mean periods therefrom. The clus- 
ters that were included in their survey turned out to 
be NGC 1835, NGC 1898, NGC 1916, NGC 1928, 
NGC 2005, and NGC 2019. For NGC 1466, we used 
the value from Walker (1992b). For the remaining 
clusters— namely, Reticulum, NGC 1786, NGC 1841, 
NGC 2210, and NGC 2257— we utilized the values 
summarized in Table 6 of Walker (1992a). Metallic- 
ity values for all LMC clusters were taken from the 
recent Mackey & Gilmore (2004b) compilation. 

An interesting characteristic found in the OGLE 
data for NGC 1835 is the presence of a bimodal pe- 
riod distribution for the RRab variables, as if the 



ftp : //ftp . astrouw. edu.pl/ogle/ ogle2/var _ stars/lmc/rrlyr/ 
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Table 2 Galactic Globular Clusters with at Least 5 Known RRab Variables 



Name Other [Fe/H] (P ab ) iV ab Type C Population a 

(d) 



NGC 362 




-1. 


.16 





.564 


13 


Ool 


-0. 


87 


YH 


NGC 1261 




-1. 


.35 





.555 


13 


Ool 


-0. 


71 


YH 


NGC 2419 




-2. 


12 





.655 


24 


OoII 


+0. 


86 


OH 


NGC 3201 




-1. 


58 





.554 


72 


Ool 


+0. 


.08 


YH 


NGC 4590 


M68 


-2. 


06 





.622 


14 


OoII 


+0. 


17 


YH 


NGC 4833 




-1. 


.88 





.708 


7 


OoII 


+0. 


.93 


OH 


NGC 5024 


M53 


-1. 


.99 





.649 


29 


OoII 


+0. 


81 


OH 


NGC 5053 




-2. 


.29 





.672 


5 


OoII 


+0. 


52 


YH 


NGC 5139 


oj Cen 


-1. 


62 





.651 


76 


OoII 


+0. 


89 


ujC 


NGC 5272 


M3 


-1. 


.57 





.555 


145 


Ool 


+0. 


18 


YH 


NGC 5466 




-2. 


.22 





.646 


13 


OoII 


+0. 


58 


YH 


NGC 5824 




-1. 


.85 





.624 


7 


OoII 


+0. 


79 


OH 


NGC 5904 


M5 


-1. 


.27 





.551 


91 


Ool 


+0. 


31 


OH 


NGC 5986 




-1. 


.58 





.652 


7 


OoII 


+0. 


97 


OH 


NGC 6121 


M4 


-1. 


.20 





.553 


31 


Ool 


-0. 


06 


OH 


NGC 6171 


M107 


-1. 


.04 





.538 


15 


Ool 


-0. 


.73 


OH 


NGC 6229 




-1. 


43 





.553 


30 


Ool 


+0 


24 


YH 


NGC 6266 


M62 


-1. 


.29 





.548 


131 


Ool 


+0. 


55 


OH 


NGC 6284 




-1. 


.32 





.588 


6 


Oo-Int 


+0. 


88 


OH 


NGC 6333 


M9 


-1. 


.75 





.638 


8 


OoII 


+0. 


87 


OH 


NGC 6341 


M92 


-2. 


28 





.630 


11 


OoII 


+0. 


.91 


OH 


NGC 6362 




-0. 


95 





.547 


18 


Ool 


-0. 


58 


OH 


NGC 6388 




-0. 


60 





.676 


9 


OoIII 


-0. 


.69 


BD 


NGC 6402 


M14 


-1. 


.39 





.564 


39 


Ool 


+0. 


65 


OH 


NGC 6426 




-2. 


.26 





.704 


9 


OoII 


+0. 


58 


YH 


NGC 6441 




-0. 


.53 





.756 


43 


OoIII 


-0. 


73 


BD 


NGC 6558 




-1. 


.44 





.556 


6 


Ool 


+0. 


70 


OH 


NGC 6584 




-1. 


49 





.560 


34 


Ool 


-0. 


.15 


YH 


NGC 6626 


M28 


-1. 


45 





.577 


8 


Ool 


+0. 


.90 


OH 


NGC 6642 




-1. 


.35 





.544 


10 


Ool 


-0. 


04 


YH 


NGC 6656 


M22 


-1. 


64 





.632 


10 


OoII 


+0. 


.91 


OH 


NGC 6712 




-1. 


.01 





.557 


7 


Ool 


-0. 


62 


OH 


NGC 6723 




-1. 


12 





.541 


23 


Ool 


-0. 


.08 


OH 


NGC 6864 


M75 


-1. 


16 





.587 


25 


Oo-Int 


-0. 


07 


OH 


NGC 6934 




-1. 


54 





.574 


68 


Ool 


+0 


25 


YH 


NGC 6981 


M72 


-1. 


.40 





.547 


24 


Ool 


+0. 


.14 


YH 


NGC 7006 




-1. 


.63 





.569 


53 


Ool 


-0. 


.28 


YH 


NGC 7078 


M15 


-2. 


.26 





.637 


39 


OoII 


+0. 


.67 


YH 


NGC 7089 


M2 


-1. 


62 





.725 


17 


OoII 


+0. 


92 


OH 


IC 4499 




-1. 


60 





.580 


63 


Oo-Int 


+0. 


.11 


YH 


Ruprecht 106 




-1. 


.67 





.617 


13 


Oo-Int 


-0. 


82 


YH 



"YH = "Young Halo"; OH = "Old Halo"; BD = "Bulge/Disk"; ujC = u> Ccntauri. From Mackey & van den Bergh (2005). 
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Table 3 Milky Way Dwarf Galaxy Satellite Globular Clusters with at Least 5 RRab Variables 



Name 


Other [Fe/H] 




iV ab 


Type 


C 


Population 






(d) 










LMC Globular Clusters 


NGC 1466 


-2.17 


0.581 


19 


Oo-Int 


+0.42 


LMC 


NGC 1786 


-1.87 


0.677 


17 


OoII 


+0.39 


LMC 


NGC 1835 


-1.79 


0.606 


55 


Oo-Int 


+0.57 


LMC 


NGC 1841 


-2.11 


0.676 


18 


OoII 


+0.71 


LMC 


NGC 1898 


-1.37 


0.565 


17 


Ool 


+0.03 


LMC 


NGC 1916 


-2.08 


0.729 


6 


OoII 


+0.97 


LMC 


NGC 1928 


-1.27 


0.634 


6 


OoII 


+0.94 


LMC 


NGC 2005 


-1.92 


0.668 


6 


OoII 


+0.90 


LMC 


NGC 2019 


-1.81 


0.606 


24 


Oo-Int 


+0.66 


LMC 


NGC 2210 


-1.97 


0.598 


13 


Oo-Int 


+0.65 


LMC 


NGC 2257 


-1.63 


0.578 


17 


Ool 


+0.42 


LMC 


Reticulum 


-1.66 


0.559 


16 


Ool 


+0.00 


LMC 


Fornax dSph Globular Clusters 


Fornax GC1 


-2.05 


0.611 


5 


Oo-Int 


-0.30 


Fornax 


Fornax GC2 


-1.83 


0.574 


15 


Ool 


+0.50 


Fornax 


Fornax GC3 


NGC 1049 -2.04 


0.606 


13 


Oo-Int 


+0.44 


Fornax 


Fornax GC4 


-1.90 


0.600 


11 


Oo-Int 


-0.42 


Fornax 


Fornax GC5 


-1.90 


0.576 


7 


Ool 


+0.52 


Fornax 


Sagittarius dSph (Candidate) Globular Clusters 


NGC 4147 


-1.83 


0.525 


12 


Ool 


+0.55 


Sag? 


NGC 5634 


-1.88 


0.660 


12 


OoII/Int 


+0.91 


Sag/OH? 


NGC 6715 


M54 -1.58 


0.590 


55 


Oo-Int 


+0.54 


Sag 


Arp 2 


-1.76 


0.584 


8 


Oo-Int 


+0.53 


Sag 




CMa dSph (Candidate) Globul 


%r Clusters 






NGC 1851 


-1.22 


0.571 


21 


Ool 


-0.32 


CMa/OH? 


NGC 2808 


-1.15 


0.563 


10 


Ool 


-0.49 


CMa/OH? 


NGC 5286 


-1.67 


0.630 


29 


OoII 


+0.80 


CMa/OH? 
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cluster actually presented an internal Oosterlioff di- 
chotomy (see Fig. 4 in Soszynski et al. 2003). A sim- 
ilar phenomenon has recently been discovered in the 
Fornax dSph globular cluster 4 (Greco et al. 2005). 
There are preliminary indications (Vidal et al. 2009, 
in preparation) that NGC 6284 may, in spite of the 
cluster's Oosterhoff-intermediate classification indi- 
cated by (Pah), lack RR Lyrae stars in the Oosterlioff 
gap region as well. 

• In the special case of the OoIII globular clusters 
NGC 6388 and NGC 6441, the newly revised (P ab ) 
values from Corwin et al. (2006) were adopted. 

• For the dSph galaxies, the values come from Table 1 
in Catclan (2004b). 

• For the SMC, we adopted the (P ab ) value from 
Soszynski et al. (2002), and a metallicity that is an 
average between the Smith et al. (1992) and Butler, 
Demarque, & Smith (1982) values. 

• For the LMC, the adopted RR Lyrae metallicity rep- 
resents an average between the spectroscopic mea- 
surements of Borissova et al. (2004) and Gratton 
et al. (2004a). The question regarding the LMC's 
(-Pab) value is more complicated: as it happens, the 
MACHO (Alcock et al. 1996) and OGLE (Soszynski 
et al. 2003) surveys provide results that differ by 
0.01 d from one another, the MACHO result, (P ab ) = 
0.583 d, placing the galaxy inside the Oosterhoff gap 
region in Figure 5, but the OGLE result, (P a b) = 
0.573 d, indicating an Ool classification instead. 
Given the large databases used in both studies — 
5455 RRab's in the case of OGLE (Soszynski et al. 
2003) and 6158 in the case of MACHO (Alcock et al. 
2003) — this difference appears to be intrinsic. This 
could be due to population gradients in the LMC, 
since whereas the OGLE fields are concentrated along 
the bar of the LMC, the MACHO survey includes 
many LMC halo fields. 13 A similar segregation has 
been suggested to exist in our own galaxy, in the 
sense that OoII (i.e., longer-period) variables would 
be more confined to relatively small distances from 
the Galactic plane, whereas Ool variables would lie 
farther from the plane on average 14 — only that, in the 



13 For a comparison between the two teams' 
area coverage, see the maps available at 

http: / /bulge .princeton. edu/ ^ogle/ogle2/rrlyr_ lmc_ map. html (OGLE 
team) and http : //wwwmacho .mcmaster . ca/Systems/Coords/LMC _ Fields . gif 

(MACHO project). 

14 This presumably implies that the HB morphology of the field 
is bluer near the plane than farther away, which is consistent 
with the results from Preston, Schcctman, & Beers (1991) and 
Kinman, Suntzcff, & Kraft (1994). Interestingly, Lee & Carney 
(1999a) note that the results of Layden (1996) for RR Lyrae 
stars and Wilhelm et al. (1996) for HB stars in general support 
an accompanying change in kinematic behavior, from progradc 



case of the LMC, the longer-period variables would 
be located preferentially farther away from the bar. 
For our present purposes, we decided to include two 
separate datapoints for the LMC RR Lyrae in Fig- 
ure 5, to reflect the possible presence of two different 
populations. We note, in addition, that Alcock et 
al. (1996) already called attention to the Oosterhoff- 
intermediate nature of the RR Lyrae period distribu- 
tion in that galaxy. For a recent discussion of the ev- 
idence for a genuine halo component in the LMC (as 
indicated primarily by the velocity dispersion of the 
RR Lyrae), see Minniti et al. (2003) and Borissova 
et al. (2004), but also Gallart et al. (2004), Zarit- 
sky (2004), and Carrera et al. (2008) (and references 
therein) for possible alternative interpretations. 

Figure 6 (right panel) compares the (P a b) distribu- 
tions for Galactic globular clusters and nearby extra- 
galactic systems. Simple eye inspection clearly reveals 
that the two distributions arc remarkably different, 
with the one for extragalactic systems strongly peaking 
at a (P a b) — 0.6 d, which however corresponds to the 
minimum of the Galactic globular cluster distribution. 
A Kolmogorov-Smirnov test confirms that the differ- 
ence is highly significant, with a Pks = 0.3338, imply- 
ing a probability that the two sets are derived from the 
same parent distribution of only Pks = 1.8%. Remov- 
ing the dwarf galaxies (but not their associated globular 
clusters) from the sample changes these figures slightly; 
we now find Pks = 0.2856, implying a Pks = 12.5%. 
(These figures do not change significantly if we keep 
only the "young halo" component among the Galactic 
globulars.) Removing NGC 6388 and NGC 6441 from 
the Galactic globular cluster sample so that the com- 
parison is limited to globular cluster systems having 
comparable metallicities, we find Pks = 0.3034, imply- 
ing a Pks = 9.2%. Finally, by assuming that ui Cen, 
NGC 1851, NGC 2808, and NGC 5286 are all bona- 
fide members of the Galactic globular cluster system — 
which is equivalent to assuming that the CMa dSph 
does not possess its own globular cluster system, which 
might perhaps be supported by the fact that the sug- 
gested CMa galaxy appears to be much too young and 
metal-rich (e.g., Sbordone et al. 2005b; Carraro, Moit- 
inho, & Vazquez 2008) to host several old and metal- 
poor globulars — we find Pks = 0.3714, with a corre- 
sponding Pks = 2.8%. 

In summary, it appears quite clear that the distribu- 
tion of mean RRab periods for the Galactic and nearby 
extragalactic systems have only a very small probability 
of being compatible with one another. 



systemic rotation to retrograde rotation as one moves away from 
the plane. 
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Fig. 6 (Left panel) The right panel of Figure 5 is repeated here for convenience. (Right panel) Histogram of {Pah) values for 
Galactic globular clusters (hatched bars) and for the nearby extragalactic sysems (gray bars). There is a marked difference 
between the two histograms, with the objects of extragalactic origin showing a peak at (P a b) — 0.6 d, which however 
corresponds to a minimum in the corresponding Galactic globular cluster distribution. 



6.2 On the Origin of the Oosterhoff Dichotomy 

In Figure 7, we show the variation in the C parame- 
ter with metallicity for Galactic globular clusters. To 
produce this figure, the main source of HB morphol- 
ogy parameters was the compilation by Mackey & van 
den Bergh (2005) for Galactic globulars, and the one 
by Mackey & Gilmore (2004b) for nearby extralactic 
globular clusters. The [Fe/H] values came from the 
Feb. 2003 edition of the Harris (1996) catalog. We 
note, however, that some of the values provided in the 
Mackey & van den Bergh compilation appear to be in- 
correct. For NGC 6388 and NGC 6441, in particular, 
they quote values of C = —1.0, which cannot be cor- 
rect given the well-known presence of extended blue HB 
components (Piotto et al. 1997; Rich et al. 1997) and 
RR Lyrae variables (Silbermann et al. 1994; Layden et 
al. 1999; Pritzl et al. 2000, 2001, 2002b, 2003; Corwin et 
al. 2006) in both clusters. Accordingly, we computed 
new C values for these globulars, from the photome- 
try provided in the quoted papers and in Piotto et al. 
(2002). Number counts for non- variable stars in the Pi- 
otto et al. study were kindly provided by M. Zoccali 
(2003, priv. comm.), whereas the number counts for 
RR Lyrae variables were taken from the Pritzl et al. 
and Corwin et al. studies. We obtain C ~ —0.69 for 
NGC 6388, and C ~ -0.73 for NGC 6441. 



In the case of M3, Catelan, Ferraro, & Rood (2001b) 
and Catelan, Rood, & Ferraro (2002b) have shown that 
the HB morphology is not constant with radius, the HB 
type closer to the center being significantly bluer than 
the HB type farther away. Previous values of the £ 
parameter for the cluster tended to be based on pho- 
tographic measurements that were carried out for the 
outer parts, so that the HB type was accordingly un- 
derestimated. When properly taking into account HB 
stars from all radial regions in the cluster in the num- 
ber counts, a value C = 0.18 obtains (Catelan 2004a) — 
which is about 0.1 bluer than provided in the avail- 
able compilations (e.g., Harris 1996; Mackey & van den 
Bergh 2005). 

u Centauri (NGC 5139) lacks an HB type measure- 
ment in both the Harris (1996) and Mackey & van 
den Bergh (2005) compilations. We accordingly adopt 
the HB type provided by Borkova & Marsakov (2000), 
namely C = 0.89. 

In Figure 7 are also displayed isochrones based on 
the work of Catelan & de Freitas Pacheco (1993). These 
are practically identical to the ones presented in their 
Figure la, the difference in age with respect to the 
middle (reference) isochrone being +0.8 Gyr (upper 



Horizontal Branch Stars: Observations and Theory 



19 



0.0 



-0.5 



.-, -1.0 
^ -1.5 



-2.0 



-2.5 



■ "Bulge, disk" GCs 
• "Old halo" GCs 
O 'Young halo" GCs 
Oosterhoff gap? 




~\ — i — i — T 



"i — i — i — r 



"i — i — i — r 



"i — i — i — r 





Oosterhoff gap? 





(o Centauri 


A 


LMC GCs 


o 


Fornax dSph GCs 




Sag. dSph GCs 


• 


CMa dSph GCs? 




-1.0 -0.5 0.0 0.5 
(B-R)/(B+V+R) 



1.0 



-0.5 0.0 0.5 
(B-R)/(B+V+R) 



1.0 



Fig. 7 (Left panel) Position of the Galactic globular clusters with a defined Oosterhoff type in the metallicity-"HB type" 
plane. The symbols are the same as in Figure 5. The region marked as a triangle and termed "Oosterhoff gap?" represents 
a seemingly "forbidden region" for bona-fide Galactic globular clusters. The overplotted lines are isochrones from Catelan 
& de Freitas Pacheco (1993). (Right panel) To the previous plot the position of the globular clusters which have been 
associated with dwarf satellite galaxies of the Milky Way are added. Only globulars with determined HB types are shown. 
Filled symbols for the extragalactic systems indicate an Oosterhoff -intermediate status. Of the 9 extragalactic globulars 
classified as Oosterhoff-intermediate, 7 fall inside the triangular region, which we tentatively call "Oosterhoff gap" region in 
analogy with Figure 5. Bona-fide Galactic globular clusters with Oosterhoff-intermediate status are marked with crosses. 



isochrone) and —2.0 Gyr (lower isochrone). 15 These 
isochrones were computed assuming that the total mass 
loss on the RGB does not depend on metallicity, as re- 
cently suggested by Origlia et al. (2002). On the other 
hand, an RGB mass loss dependence on metallicity, 
as suggested by the expressions shown in Table 1 (see 



5 Note, however, that those authors also pointed out that there is 
a "degeneracy effect" associated with the ages of these isochrones, 
since several different ages can provide morphologically indistin- 
guishable isochrones, depending on the assumptions regarding 
the amount of mass loss on the RGB and the chemical composi- 
tion. Conversely, different assumptions for the chemical composi- 
tion and RGB mass loss can be adopted to mimic a favored rela- 
tive age scale on the C — [Fe /H] plane [which is actually the reason 
why this particular set of isochrones was chosen in the present 
study: it appears to provide a reasonable analog of Fig. 9b in Rey 
ct al. (2001) in the scenario in which mass loss does not depend 
on metallicity (Origlia ct al. 2002)] . Accordingly, we expect that 
different "age labels" may be assigned to these same isochrones 
in the future, reflecting (for instance) changes in the absolute age 
and/or RGB mass loss scales for Galactic globular clusters. Note 
also that the WMAP and SDSS results on the age of the Universe 
(Spergel ct al. 2003, 2007; Tegmark ct al. 2004, 2006) strongly 
constrain the range of possible solutions among those discussed 
by Catelan & de Freitas Pacheco (1993). 



Fig. 4), would lead to somewhat distorted isochrones 
in the C — [Fe/H] plane, predicting redder HB types at 
lower [Fe/H] (and vice- versa) than would be the case 
for a similar, constant-AM isochrone. This effect is 
clearly shown in Fig. 9 of Rey et al. (2001), which com- 
pares isochrones computed with constant mass loss and 
isochrones computed assuming a mass loss rate as given 
by the Reimers (1975a, b) mass loss formula. 

Note that these isochrones can be well described by 
a modified Fermi-Dirac profile as follows: 

C= -, r--l. (1) 

l +7 ex P (M±^) 

The middle isochrone provides a reasonable approxi- 
mation to the inner-halo globular clusters, which in 
fact comprise a large fraction of the "old halo" glob- 
ulars shown in Figure 5. It is also very similar to 
the isochrones labeled "At = 0.0" in Rey et al. (2001) 
(their Fig. 9). For this isochrone, the parameters used 
are a = 1.227, f3 = 0.130, and 7 = 1. For the 
lower isochrone, which in Catelan & de Freitas Pacheco 
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(1993) corresponded to the LMC globular clusters, one 
has a — 1.707, [3 — 0.140, and 7=1; the relative age 
scale is fairly similar to that in Fig. 9b of Rey ct al., 
except (as expected) for the most metal-poor clusters 
with the redder HB types. Finally, the upper isochrone 
represents the RR Lyrae stars in the Galactic bulge, 
and is well described by a = 1.027, (3 = 0.130, and 
7 = 1. 

Note that the plots in Figure 7 actually use an in- 
verted form of eq. (1), which we also provide for conve- 
nience: 



[Fc/H] = [3 In 



1 + jC 



1 



— a. 



(2) 



Figure 7 (left panel) shows very clearly that there 
is a zone of avoidance for RR Lyrae-rich globular clus- 
ters in our galaxy, defined by a rectangle located in 
the region -2.18 < [Fe/H] < -1.5, 0.3 < C < 0.75. 
However, in the right panel, when globular clusters as- 
sociated with the neighboring galaxies are plotted, one 
sees that this region is actually preferentially occupied 
by the external globulars. The situation immediately 
brings to mind the similar phenomenon that was found 
in our analysis of Figure 5; in that case, the "zone of 
avoidance" for the Galactic globulars — the Oosterhoff 
gap — was also the zone preferentially occupied by the 
globulars associated with neighboring galaxies. There- 
fore, one cannot help but suspect that the noted "avoid- 
ance region" in the left panel of Figure 7 will also be 
related to the Oosterhoff-intermediate globulars, as was 
the Oosterhoff gap region in Figure 5. Indeed, Fig- 
ure 7 (right panel) clearly shows that nearby extra- 
galactic globulars do not only preferentially fall in the 
quoted region of this diagram: there is in fact a smaller, 
well-defined, triangular-shaped region in the plot where 
most Oosterhoff-intermediate globulars of extragalactic 
origin (7 out of 9) can be found. The vertices of this 
triangle, in the [Fe/H] — £ plane, are given by the fol- 
lowing coordinates: (—1.5, 0.525); (—1.92, 0.76); and 
(-2.26, 0.36). 

Note that the four Oosterhoff-intermediate Galac- 
tic globulars (NGC 6284; M75 = NGC 6864; IC 4499; 
Ruprccht 106), indicated by large "x" symbols in the 
right-hand plot, arc located very far away from the 
"Oosterhoff gap" region in this diagram, and may ac- 
cordingly represent different phenomena, or even be 
due to statistical fluctuations (see also Catelan 2004b) . 
Note, in this sense, that NGC 6284 has a mere 6 known 
RRab's (Clement et al. 2001), none of which with a 
period that falls in the Oosterhoff gap range (i.e., 0.58- 
0.62 d); that M75 has a multimodal HB; and that the 
RRab's in Rup 106 (which entirely lacks RRc's) all fall 



very close to the red edge of the instability strip. Fi- 
nally, IC 4499 is at the very edge of the Oosterhoff gap 
region, with a mean ab-type period of exactly 0.580 d, 
and might as well have been classified as Ool. 

This scenario provides unprecedented detail about 
the origin of the Oosterhoff dichotomy, going signifi- 
cantly beyond the original and insightful early analyses 
of the problem by Castellani (1983), Rcnzini (1983), 
Lee & Zinn (1990), and Bono, Caputo, & Stcllingw- 
erf (1994) — who previously suggested that the origin of 
the Oosterhoff dichotomy among Galactic globular clus- 
ters was the lack of RR Lyrae-rich globulars (or, more 
specifically, the presence of a majority of clusters with 
exclusively blue HBs) at metallicities intermediate be- 
tween the bulk of the Ool and the OoII clusters. Indeed, 
theoretical predictions of an "Oosterhoff-intermediate 
area" in the [Fe/H] — C plane have previously been pro- 
vided by Lee & Zinn and Bono et al.; the results of 
these studies are summarized in Figure 8, which shows 
reasonable success in predicting which globulars should 
be Oosterhoff-intermediate. It is interesting to note 
that the two Oosterhoff-intermediate Fornax dSph glob- 
ular clusters which lie farther away from the triangular- 
shaped region in Figure 8, as well as Rup 106, fall very 
close to the Oosterhoff-intermediate region located to- 
wards red HB types in this plane, in agreement with 
the theoretical predictions by Lee & Zinn. (A region 
equivalent to the area labeled "Bono et al." in this plot 
was also predicted, though over a somewhat more lim- 
ited C range for a given [Fe/H], by Lee & Zinn — see 
their Fig. 2.) 

An important question that remains unanswered and 
which will certainly require additional work is the fol- 
lowing: why did the Galactic globular cluster system 
avoid the triangular-shaped region in Figure 7, whereas 
the nearby extragalactic globulars were instead prefer- 
entially "attracted" to it? According to Figure 7, part 
of the explanation could be that the Galactic "old halo" 
and nearby extragalactic globular clusters differ in age 
by ~ 2 Gyr (on average). However, little (if any) ev- 
idence for an age difference has been found so far: in 
the case of the Fornax dSph, Buonanno et al. (1998) 
show that all globular clusters have similar ages, com- 
parable to those found in Galactic globular clusters of 
similar metallicity. The exception appears to be clus- 
ter 4, which is ~ 3 Gyr younger than the other For- 
nax dSph globulars (Buonanno et al. 1999), and which 
indeed appears to have the redder HB type (see Ta- 
ble 3). In the case of the LMC, virtually all "old" 
clusters have ages comparable to those of the oldest 
Galactic globular clusters, a point which has been re- 
peatedly emphasized in the recent literature (e.g., Bro- 
cato et al. 1996; Johnson ct al. 1999, 2002; Mackey & 
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Fig. 8 Same as in the right panel of the previous figure, but 
now showing, as hatched regions, the predicted Oosterhoff- 
intermediate region, according to the model calculations by 
Bono et al. (1994, right, in light gray) and Lee & Zinn (1990, 
left, in dark gray). 



Gilmorc 2004a). Among the Oosterhoff- intermediate 
globular clusters in the Sagittarius dSph, only Arp 2 
was previously suggested to be "young" (but note that 
it has a predominantly blue HB), although the recent 
VandenBerg (2000) and Salaris & Weiss (2002) studies 
suggest that its age is probably not much lower than 
those of Galactic globular clusters of similar mctallic- 
ity. 16 Also, Layden & Sarajedini (2000) argue that the 
Sagittarius dSph globulars M54 and Arp 2 (as well as 
Ter 8) are basically coeval with the bulk of the Galactic 
halo globular clusters. Therefore, it appears difficult, 
given the available age determinations for Oosterhoff- 
intermediate globular clusters, to argue that they differ 
in age significantly with respect to OoII and Ool globu- 
lars. In a similar vein, it should be noted that the sepa- 
ration of the clusters into "young halo" and "old halo" 
groups by Mackey & van den Bergh (2005) is some- 
what artificial, since there is clearly significant overlap 
in turnoff ages between the two groups, as can be seen 
from Figure 9 in Mackey & Gilmore (2004b): according 



16 Among the OosterhofT-intermediatc Galactic globular clusters, 
both IC 4499 (Fcrraro et al. 1995) and Rup 106 (Buonanno ct 
al. 1993) were indeed previously found to be young globular clus- 
ters, although the case for a young IC 4499 is not supported by 
VandenBerg (2000), Salaris & Weiss (2002), or De Angeli ct al. 
(2005), Rup 106 being instead more consistent, according to these 
studies, with an age slightly lower than found for the bulk of the 
Galactic globular clusters. 



to the quoted figure, the oldest "young halo" globular 
clusters are at least as old as the oldest "old halo" clus- 
ters, and likewise, the youngest "old halo" clusters have 
ages that are similar to those of the youngest "young 
halo" clusters. This clearly shows the perils of using 
HB morphology as an age indicator for individual glob- 
ular clusters. We note, in passing, that De Angeli et 
al. (2005) find no clear correlation between the ages of 
Galactic globular clusters and their galactocentric dis- 
tances. 

6.3 Implications for the Formation History of the 
Galaxy 

Several authors have suggested that the "young halo" 
Galactic globular cluster system at least may have been 
formed from the accretion of dwarf galaxies resembling 
the present-day dwarf satellites of the Milky Way (e.g., 
Searle & Zinn 1978; Zinn 1993b; Mackey & van den 
Bergh 2005). However, we have already seen (Sect. 6.1) 
that the distribution of RR Lyrae periods for the Galac- 
tic globulars, including the "young halo" component, is 
different from that of those systems. It is instructive to 
perform a similar quantitative comparison between the 
distribution of HB types for the Galactic and nearby 
extragalactic globular clusters. 

Using only the "young halo" component, we find a 
Kolmogorov-Smirnov value of Dks = 0.5425, implying 
a probability of only Pks = 0.2% that the "young halo" 
Galactic globulars and the nearby extragalactic globu- 
lars were drawn from the same parent population. If we 
remove the LMC globulars from the sample (given that 
many authors suggest that the protogalactic fragments 
that led to the formation of the Milky Way were ac- 
tually dSph-like, not LMC-like), we find D KS = 0.498, 
implying a probability of only Pks = 2.4% that the 
samples were drawn from the same parent distribution. 
Therefore, the Kolmogorov-Smirnov test rejects, with 
high statistical significance, the hypothesis that the 
Galactic globular cluster system was assembled from 
the capture of protogalactic fragments that resembled 
the present-day dwarf satellite galaxies of the Milky 
Way — Fornax and Sagittarius in particular. 

This does not mean, of course, that there may not 
have been a (now extinct) primordial dwarf galaxy pop- 
ulation with properties that were different from those of 
the surviving dwarf satellite galaxies of the Milky Way. 
We thus face the task of defining whether our current 
sample of dwarf galaxies is atypical in some respect. 
For instance, it is now well established that the detailed 
chemical patterns of even the more metal-poor popu- 
lations of the Milky Way dwarf satellite galaxies bear 
little resemblance to that found among most stars in 
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our galaxy (e.g., Venn et al. 2004). Note, however, that 
our satellite galaxies do not appear to be anomalous in 
this respect, since Bonifacio et al. (2004; see especially 
their Fig. 5) have shown that their abundance patterns 
are consistent with those for objects which are sup- 
posed to have given rise to dSph galaxies — namely, the 
Damped Lyman a (DLA) systems (e.g., Matteucci, Mo- 
laro, & Vladilo 1997; Hopkins, Rao, & Turnshck 2005). 
Also, the globular clusters which were unquestionably 
accreted by the Milky Way from dSph fragments and 
which have been carefully studied for chemical abun- 
dance patterns, such as M54 (Brown, Wallcrstcin, & 
Gonzalez 1999), Palomar 12 (Brown et al. 1997; Cohen 
2004), and Terzan 7 (Tautvaisicnc et al. 2004; Sbor- 
done et al. 2005a, 2007), show clearly different abun- 
dance patterns with respect to bona-fide Galactic glob- 
ulars. A careful and detailed discussion of this point 
has recently been presented by Pritzl, Venn, & Irwin 
(2005b). In fact, the presence of such peculiar abun- 
dance patterns, particularly a lowered abundance of the 
alpha-captured elements, has recently been viewed as a 
signature of a possible accretion origin; the reader is 
referred to the recent spectroscopic study of NGC 5694 
by Lee, Lopez-Morales, & Carney (2006) for a recent 
example. Note, on the other hand, that there are at 
least some outer-halo globular clusters with completely 
normal abundance patterns (Ivans et al. 2001; Cohen 
& Melendez 2005; see also footnote 19 below), thus sug- 
gesting that it is not only the inner Galactic halo that 
contains globular clusters that formed "in situ." On the 
other hand, Mottini, Wallcrstcin, & McWilliam (2008) 
have recently claimed that Arp 2 and Terzan 8 — which 
are also Sagittarius globulars — both have alpha-element 
abundances that are typical of Galactic halo clusters, 
which certainly confuses the chemical picture some- 
what. As far as the newly discovered SDSS dSph galax- 
ies (e.g., Belokurov et al. 2006, 2007), recent surveys 
have revealed that at least some of them are also Oost- 
crhoff intermediate (e.g., Kuehn et al. 2008) and present 
unusual abundance patterns (Frebcl et al. 2009, in 
preparation), although seemingly "regular" OoII dwarfs 
are also present (Siegel 2006; Dall'Ora et al. 2006; 
Greco et al. 2008; Musella et al. 2009, in preparation; 
see also Catelan 2009 for a critical discussion). 

In conclusion, most of the recent evidence suggests 
that, even though we seem to be witnessing mergers 
between dSph satellite galaxies and our own galaxy "in 
real time" (e.g., Sagittarius, CMa), these are not truly 
representative of the events that led to the formation of 
the Milky Way, and only a relatively minor fraction of 
the Galactic halo may have been assembled from dSph- 
likc protogalactic fragments resembling the present-day 
Milky Way dSph satellites (see also Geisler et al. 2007 
for a recent discussion). 



6. 3. 1 RR Lyrae Stars in MSI and the Origin of the 
Galactic Halo 

It has recently been suggested that the "young" second- 
parameter globular clusters in our galaxy (as well as 
at least some of the Milky Way dSph satellites) were 
accreted from M31 when the latter was forming its 
Population II stars (Kravtsov 2002). Given that our 
galaxy's globular clusters clearly present the Ooster- 
hoff dichotomy, but not so its dSph satellite galaxies, 
the question naturally presents itself: how does M31 
classify, in terms of Oosterhoff status? 

Of course, this is a very difficult question to answer, 
given that the detection of RR Lyrae variable stars at 
the distance of M31 is very difficult from the ground. It 
was originally attempted by Pritchet & van den Bcrgh 
(1987) using the CFH 3.6m telescope. These authors 
claimed the detection of 30 probable RR Lyrae stars in 
their surveyed field, thus suggesting that the M31 halo 
has a specific frequency of RR Lyrae variables compara- 
ble to that in RR Lyrae-rich Galactic globular clusters. 
They also reported a mean period for the ab-type vari- 
ables of (P a b) = 0.548 d, thus classifying the M31 halo 
field as Ool. 

More recently, Dolphin et al. (2004) performed a 
variability survey, using the WIYN 3.5m telescope, of 
an M31 halo field that includes the Pritchet & van den 
Bcrgh (1987) field, finding a specific frequency of RR 
Lyrae variables dramatically lower than in the Pritchet 
& van den Bergh study. Using time-series observations 
with ACS onboard the HST, Brown et al. (2004a) have 
shown that, in fact, the specific RR Lyrae frequency of 
the M31 halo is intermediate between the two quoted 
studies. Interestingly, they also find a (P a b) = 0.594 d, 
which corresponds to the Oosterhoff-intermediate do- 
main in Figure 5. On the other hand, as we have seen, 
the dSph satellites of the Milky Way are indeed almost 
exclusively Oosterhoff-intermediate — and so are prob- 
ably the dSph satellites of M31 as well (Pritzl et al. 
2002a, 2005a). Therefore, the incorporation of dSph 
satellites from the M31 system into our own galaxy's 
would not fundamentally change the stellar pulsation 
properties of the Milky Way's dSph system. 

To be sure, extensive variability surveys of not only 
the M31 halo field and its dwarf satellite galaxies, but 
also (and most challenging of all) of its globular clus- 
ters, are badly needed to place constraints on its old 
halo's (and indeed, as we have just seen, the whole Lo- 
cal Group's) formation history. It would be very im- 
portant to establish whether the Oosterhoff dichotomy 
is present in the M31 globulars, as well as in M31's 
general field, since this would enable a systematic com- 
parison between the very oldest stars in the M31 dSph 
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satellites and their peers in the general M31 halo, thus 
posing constraints on the extent to which the latter 
may have been assembled from "building blocks" sim- 
ilar to the former. Unfortunately, such studies have 
so far remained rather limited (Clemcntini ct al. 2001; 
Contreras et al. 2008), given that they require expensive 
time-series observations from space. Ground-based sur- 
veys could in principle be carried out using large tele- 
scopes equipped with adaptive optics instruments, but 
unfortunately, this technique remains limited to near- 
infrared observations — a wavelength regime where not 
only the red giants are brightest but also the ampli- 
tudes of RR Lyrae stars are smallest (e.g., Longmore et 
al. 1985; Liu & Janes 1990b; Jones, Carney, & Fulbright 
1996), thus making their detection hardest. 

7 The Second-Parameter Problem 

The second parameter problem of globular cluster as- 
tronomy is commonly defined as the existence of glob- 
ular clusters with very similar metallicity (the "first 
parameter"; Sandage & Wallerstein 1960) and yet dif- 
ferent HB morphologies. It appears to have been first 
noted by Sandage & Wallerstein (see pages 607 and 608 
in their paper), and later also by Sandage & Wildey 
(1967) and van den Bergh (1967). The phenomenon 
has now been recognized to exist also in the M31 halo 
(Mackey et al. 2006), as well as in the Fornax dSph 
galaxy (Buonanno et al. 1998, 1999). 

The second parameter candidate that was first noted 
in the literature seems to have been age. We quote 
Sandage & Wallerstein (1960): 

... the character of the horizontal branch is spoiled by 
the two clusters M13 and M22. [Individual stars in] M13 ap- 
pear to be metal-rich, whereas the character of the horizon- 
tal branch simulates that of the very weak-lined group (M15, 
M92, NGC 5897). [...] M13 is younger than M2 or M5 (Arp 
1959). Consequently, in addition to chemical composition, the 
second parameter of age may be affecting the correlations... 

(Note that the sense of the claimed age difference be- 
tween M13 and other clusters of similar metallicity but 
redder HB type is the opposite of what modern stellar 
evolution theory indicates to be necessary to account 
for their differences in HB types.) 

The next candidate second parameter was the helium 
abundance (Sandage & Wildey 1967; van den Bergh 
1967). Later on, with the advent of the Searle & 
Zinn (1978) scenario for the formation of the Galactic 
halo — which was based upon the hypothesis that age 
is the second parameter, outer-halo globular clusters 
with predominantly red HBs being younger than inner- 
halo globulars with bluer HBs (see their Fig. 10) — 
age has quickly become the most popular second pa- 
rameter candidate. However, the age hypothesis has 



proven controversial, with different perspectives hav- 
ing been presented on both the presence of significant 
age differences among Galactic globular clusters (e.g., 
Stetson et al. 1996; Sarajedini et al. 1997; and refer- 
ences therein) and the amount that may be required 
to account for different second-parameter pairs entirely 
in terms of age (e.g., Catelan & de Freitas Pacheco 
1993; Lee et al. 1994; Ferraro et al. 1997). In ad- 
dition, other candidate second parameters abound in 
the literature; these include (or are at least related 
to, and often involve combinations of) mass loss (Pe- 
terson 1982; Catelan 2000), cluster ellipticity (Norris 
1983, 1987), stellar rotation (Mengel & Gross 1976; 
Fusi Pecci & Rcnzini 1978; Peterson 1985), magnetic 
fields (Rood & Scitzcr 1981; Castellani 1983; Castel- 
lani & Paterno 1984), the cyanogen distribution among 
red giants (Norris 1981; Norris ct al. 1981; Smith & 
Norris 1983), [CNO/Fe] (Rood & Seitzer 1981), supcr- 
oxygen-poor RGB stars (Catelan & de Freitas Pacheco 
1995), the Na-0 correlation among cluster giants (Car- 
retta & Gratton 1996; Carretta et al. 2007), the helium- 
core mass at the helium flash (Dcmarque, Mengel, & 
Sweigart 1972; Sweigart & Catelan 1998), helium mix- 
ing (VandcnBerg & Smith 1988; Langer & Hoffman 
1995; Sweigart 1997a,b; Sweigart & Catelan 1998; Cav- 
allo & Nagar 2000; Aikawa, Fujimoto, & Kato 2001; 
Caloi 2001), planetary systems (Soker 1998; Siess & 
Livio 1999; Soker & Harpaz 2000; Soker & Hadar 2001; 
Soker, Rappaport, & Fregeau 2001; Livio & Soker 2002; 
Soker & Harpaz 2007), globular cluster core density or 
concentration (Fusi Pecci et al. 1993; Buonanno et al. 
1997), and cluster mass (Recio-Blanco et al. 2006). A 
very informative and relatively recent review of the sev- 
eral "second parameters" that have been proposed in 
the literature is provided by Fusi Pecci & Bcllazzini 
(1997). 

Recently, Cavallo, Suntzeff, & Pilachowski (2004) 
have argued that neither the deep mixing nor the pri- 
mordial channel are capable of satisfactorily accounting 
for the abundance patterns observed among globular 
cluster red giant stars. In the same vein, Sneden et 
al. (2004) present an extremely interesting conundrum 
in regard to the extreme abundance anomalies seen in 
stars close to the RGB tip in M13, and whose solution 
may hold the key to the extent to which deep mixing 
and/or primordial contamination may affect the evolu- 
tion of the stars on the HB: 

The correlation of O and the Mg isotopic abundance ra- 
tio suggests that M13 giants in a very small interval of initial 
mass preferentially underwent severe pollution or ablation or, 
alternatively, were preferentially formed from material secu- 
larly ejected earlier from more massive cluster giants. These 
stars only just at this moment arrived essentially at the red 
giant tip, exactly where the effect of deep mixing, if it in- 
deed exists, would likely be most easily manifest. It is a most 
remarkable coincidence, if indeed that's what it is. 
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Fig. 9 Effect of an increased helium abundance upon the 
gravities of HB stars. The lines labeled "ZAHB" correspond 
to the zero-age HB, whereas the lines labeled "TAHB" cor- 
respond to the terminal-age HB (i.e., to core helium ex- 
haustion). The vertical line indicates the temperature of 
the Grundahl et al. (1999) jump. 

7.1 Helium as a Second Parameter: Some Empirical 
Constraints 

As mentioned previously, the helium abundance was 
one of the first second parameter candidates to be pro- 
posed in the literature. Recently, with observations of 
faint stars in globular clusters having shown that many 
of the abundance patterns observed among bright gi- 
ants that might have been ascribed to deep mixing 
are also present down to the main-sequence level (see 
Gratton et al. 2004b for an extensive review and refer- 
ences), thus severely constraining the extent to which 
deep mixing may contribute to the observed patterns 
among the brighter giants, the helium abundance sce- 
nario has regained impetus, now under the label of 
primordial contamination by a previous generation of 
(likely massive AGB) stars (e.g., D'Antona et al. 2002, 
2005; D'Antona & Caloi 2004, 2008; Norris 2004; Caloi 
& D'Antona 2005, 2007; Lee et al. 2005b; Busso et al. 
2007; Ncwsham & Terndrup 2007; Piotto et al. 2007; 
Yoon et al. 2008; and references therein). For critical 
discussions of different theoretical scenarios for the for- 
mation of helium-enhanced (sub-)populations in globu- 
lar clusters, the reader is referred to the recent papers 
by Bekki & Norris (2006), Chuzhoy (2006), Maeder & 
Meynet (2006), Prantzos & Charbonncl (2006), Bekki 
et al. (2007), D'Antona & Ventura (2007), Choi & Yi 
(2007, 2008), and Decressin, Charbonnel, & Meynet 
(2007). The earlier paper by Shi (1995, see its §4), 
which unfortunately has received little attention, al- 
ready contained an interesting discussion of several dif- 



ferent scenarios that could generate helium-enhanced 
stars in globular clusters. In turn, modern views of 
the deep helium mixing scenario have recently been ad- 
vanced as well, including papers by Suda & Fujimoto 
(2006), Suda et al. (2007), and Yamada, Okazaki, & 
Fujimoto (2008). These authors propose that, in the 
dense environments of globular clusters, such deep mix- 
ing can indeed take place along the RGB, as triggered 
by tidal interactions. According to them, these phe- 
nomena could explain not only some level of He en- 
richment on the HB, but also the high rotation veloci- 
ties that are observed among sufficiently cool blue HB 
stars. The same authors also suggest that cluster main 
sequence stars can likewise get polluted, by encounters 
with AGB stars (Tsujimoto, Shigeyama, & Suda 2007). 
Soker & Harpaz (2007) also speculate that He mixing 
can take place during the pre-ZAHB phase, triggered 
by the engulfment of low-mass companions and the re- 
sulting spin-up of the star's envelope. 

What are the constraints that may be posed on the 
enhanced helium hypothesis, using HB stars? In Fig- 
ure 9, we show the impact of an enhanced helium abun- 
dance upon the logg — logT e ff diagram for blue HB 
stars. These results are based on the evolutionary cal- 
culations by Sweigart & Catelan (1998). As one can 
clearly see, a large increase in the helium abundance is 
expected to lead to a marked signature in this plane, 
with lower gravities resulting due to the increased lu- 
minosities that are brought about, with an increase in 
Y, for stars with efficient H-burning shells. For stars 
cooler than the Grundahl et al. (1999) "jump" (at about 
11,500 K; see Sect. 2.4), comparison between the spec- 
troscopically derived gravities and the theoretical pre- 
dictions should be fairly straightfoward, the only note- 
worthy complication being associated with the distor- 
tion of the Balmer lines that may be effected by the pre- 
dicted presence of a stellar wind (Vink & Cassisi 2002). 
In the case of M13, an increase in Y up to 0.28 has 
been suggested (Caloi & D'Antona 2005). The spec- 
troscopic results by Moehler et al. (2003) do not rule 
out a Yms — 0.28 for these cooler stars (but note again 
that proper inclusion of stellar winds might bring the 
derived gravities up, and therefore the required Y val- 
ues down), but appear inconsistent with a Yms ~ 0.33 
or higher (see their Fig. 8b); in addition, they do not 
favor a difference in Y between Ml 3 and M3, although 
more data would certainly be needed to firmly estab- 
lish this point. Note that the high helium abundance 
suggested for blue HB stars in NGC 2808 could also be 
constrained in this way, since Lee et al. (2005b) sug- 
gest that even the cooler blue HB stars in this clus- 
ter have helium abundances as high as Yms = 0.33, 
which should leave an obvious mark in the log g— log T e g 
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plane (Fig. 9). The helium enhancement proposed for 
the same cluster by D'Antona & Caloi (2004) is signifi- 
cantly smaller, being at the level of Yjyis = 0.27 for the 
cooler blue HB stars; the difference should be testable 
with a sufficiently large sample of stars. The same tech- 
nique can likely be used to constrain the suggestion by 
Caloi & D'Antona (2008) that even in the case of M3 
the blue HB stars may present some degree of helium 
enhancement. 

The comparison between evolutionary model predic- 
tions and spectroscopically derived gravities for tem- 
peratures higher than the Grundahl jump is a much 
more complicated affair, since the presence of radia- 
tive levitation/gravitational diffusion effects in this re- 
gion require the computation of detailed model atmo- 
spheres properly taking into account the extremely com- 
plex observed abundance patterns (Bonifacio, Castelli, 
& Hack f995; Castelli, Parthasaraty, & Hack 1997) in 
order for reliable gravities, bolometric corrections and 
color transformation to be derived. To the best of our 
knowledge, no such detailed calculations have been car- 
ried out yet. Accordingly, it remains unclear whether 
models with an overall enhancement in all metals to 
super-solar levels (which is not what the quoted abun- 
dance studies indicate), as commonly employed in the 
literature, are adequate to derive gravities of stars hot- 
ter than the Grundahl jump. 

7.2 Why Arc the Most Metal-Poor Globular Clusters 
Not the Ones with the Bluest HBs? 

The influence of the "first parameter" (metallicity) 
upon HB morphology can be summarized as the ob- 
served trend of HB type getting redder as metallic- 
ity increases (Sandage & Wallerstein I960). 17 Inter- 
estingly though, the most metal-poor globular clusters 
are not the ones with the bluest observed HBs, contrary 
to what might be expected from basic theory (see the 
isochrones overplotted on Figure 7). In fact, the trend 
of HB type becoming bluer with decreasing [Fe/H] ap- 
pears to be reversed at a [Fe/H] w —1.8 (see Fig. 7, left 
panel). The reason for such a reversal of trend is not en- 
tirely clear; while some have suggested that it could be 
due to a decrease in mass loss eficiency as metallicity 
decreases (Rood 1973), others have suggested instead 



17 Fusi Pecci et al. (1993) (see their Sect. 2.1) provide an en- 
lightening discussion of the dependence of HB temperature on 
metallicity and RGB mass loss, showing how it can become very 
difficult for a metal-rich system to produce blue HB stars, com- 
pared to metal-poor systems. Note, in addition, that the RGB 
mass increases with metallicity for a fixed age and helium abun- 
dance [see, e.g., cq. (1) in Catelan & de Freitas Pacheco 1993], 
which also favors the production of redder HB types at higher Z. 
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that a majority of the metal-poor globulars may in fact 
be somewhat younger than Galactic globular clusters of 
intermediate metallicity (Yoon & Lee 2002). According 
to the isochrones shown in Figure 7, the required age 
difference between the "old halo" globular clusters and 
the most metal-poor ones with 0.1 < C < 0.7 would 
have to be fairly large. On the other hand, the HB 
morphology-based age difference would decrease using 
mass loss formulae which imply a strong mass loss de- 
pendence on metallicity, as indeed favored by the mass 
formulae listed in Table 1 (sec Fig. 4) but contrary, 
it should be noted, to the recent results by Origlia et 
al. (2002) (see Sect. 5 above). In any case, inspection 
of recent papers in which globular cluster ages have 
been computed provides little support for an age differ- 
ence, except perhaps in the case of M68 (NGC 4590) 
(Rosenberg et al. 1999; VandenBerg 2000; Salaris & 
Weiss 2002; De Angeli et al. 2005). Interestingly, M68 
does have a redder HB type than most other Galactic 
globular clusters of similar [Fe/H] — and the possibility 
that it may be relatively young was raised early on by 
Chaboyer et al. (1996). 

7.3 A Reanalysis of Specific Second-Parameter Pairs 

In what follows, we will readdress some classical second- 
parameter pairs, including NGC 288/NGC 362 and 
M13/M3, with the goal of answering the following ques- 
tion: are the measured turnoff age differences between 
these clusters sufficient to explain the observed differ- 
ence in their HB types? We will also briefly discuss the 
case of the "young" outer-halo globular clusters Pal 3, 
Pal 4, and Eridanus, and provide a summary of recent 
developments on the second-parameter effect in metal- 
rich globular clusters. Note that an RGB mass loss that 
presents no dependence on L, g, or R, as suggested by 
Origlia et al. (2002), would require larger age differ- 
ences than computed in the subsections below in order 
to account for any given second-parameter pair in terms 
of age. 

7.3.1 The Pair NGC 288-NGC 362 

This is probably the second-parameter pair which has 
attracted the most attention in the literature, given 
the very similar mctallicities of the two clusters but 
their strikingly different HB types — NGC 362 with a 
predominantly red HB, and NGC 288 with an almost 
entirely blue HB. Reviews and extensive references to 
early work have been recently provided by Bellazzini et 
al. (2001) and Catelan et al. (2001a). Here we present 
new model calculations computed along the same veins 
as in Catelan et al. , using the several different prescrip- 
tions for mass loss on the RGB summarized in Table 1 
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above, and compare the results with the age difference 
for this pair, as carefully derived by Bellazzini et al. 
on the basis of the so-called "bridge method" (Stetson 
et al. 1996). The essence of this method is that the 
HB of NGC 1851 looks very similar to the sum of the 
NGC 288 and NGC 362 HBs, so that, by carefully su- 
perposing the latter CMDs with that for NGC 1851, 
the relative positions of their turnoffs — and hence their 
relative ages — can be derived. 18 

Figure 10 shows the result of this comparison. In 
this plot, the lines indicate the age difference that is 
required, from the standpoint of canonical stellar evo- 
lution theory, to account for the observed difference in 
HB morphology between the clusters, under the as- 
sumption that the second parameter is age, and for 
the different indicated recipes for the RGB mass loss. 
The gray bands illustrate, in turn, the measured age 
difference. (These bands are not horizontal because 
the age difference depends on the absolute age value 
itself, the log of the age difference remaining instead 
more approximately constant as a function of age.) 
The upper plot corresponds to an adopted metallic- 
ity Z = 0.001, whereas the bottom plot to a metal- 
licity Z = 0.002. This figure is similar to Figure 7 
in Catelan et al. (2001a), but one sees that the new 
results require slightly larger age differences between 
NGC 288 and NGC 362 to account for their relative 
HB types. The bottom plot clearly shows that, if the 
metallicity of these clusters is of order Z — 0.002 and 
NGC 288 is younger than about 12 Gyr, the pair can 
be accounted for in terms of age as the second param- 
eter, irrespective of the mass loss formula used. Such 
a metallicity corresponds to an [Fe/H] ~ —1.2 for an 
[a/Fe] = +0.3 (Salaris, Chicffi, & Straniero 1993). Such 
an [Fe/H] value is very similar to the values provided 
in the Feb. 2003 edition of the Harris (1996) catalog, 
and only slightly lower than obtained in the Carretta 
& Gratton (1997) scale— namely, [Fc/H] = -1.07 for 
NGC 288, and [Fe/H] = -1.15 for NGC 362. There- 
fore, it appears fair to say that, in the Carretta & Grat- 
ton scale, this pair is indeed consistent with age as the 
sole second parameter. 



18 Underlying this method is of course the assumption that 
NGC 1851 is chemically similar to both NGC 288 and NGC 362, 
as discussed in detail by Bellazzini et al. (2001). The reader 
should keep in mind the possibility that this may not be strictly 
true, given the abundance peculiarities recently detected in 
NGC 1851 by Yong & Grundahl (2008), the discovery of a bi- 
modal subgiant branch in the cluster by Milonc ct al. (2008), 
and the results of the theoretical analyses by Cassisi ct al. (2008) 
and Salaris, Cassisi, & Pictrinferni (2008). Note, on the other 
hand, that fairly tight constraints on the possibility of a spread 
in Y in the cluster were recently provided by Catelan (2009). 



On the other hand, the case Z = 0.001, which cor- 
responds to an [Fe/H] ~ —1.5 for an [a/Fe] = +0.3, 
does not provide equally satisfactory consistency with 
the hypothesis that age is the second parameter. As 
can be seen from the plot, only for NGC 288 ages lower 
than about 9 Gyr is it that consistency with the rela- 
tive turnoff age difference is recovered, irrespective of 
the mass loss formulation used. For ages higher than 
10 Gyr, several mass loss formulae lead to results that 
allow one to question the hypothesis that age is the 
(only) second parameter for this pair. The best con- 
sistency with the age hypothesis obtains when using 
the "Judge-Stencel" formula; the worst, when using 
the "Goldberg" formula (see Table 1). Note that an 
[Fc/H] ~ —1.5 falls within the uncertainty range of the 
Zinn & West (1984) metallicity value, since these au- 
thors give [Fe/H] = -1.40 + 0.12 for NGC 288, and 
[Fe/H] = -1.27 ± 0.07 for NGC 362; similar values are 
provided by Kraft & Ivans (2003) and Rutledge et al. 
(1997) (in the Zinn & West scale). Therefore, we con- 
clude that it is more difficult to account for the pair 
NGC 288/NGC 362 entirely in terms of age if the Zinn 
& West scale better describes the actual abundances of 
globular cluster stars. 

It should be noted, in this sense, that the recent re- 
sults by Asplund et al. (2004) (see also Melendez 2004) 
for the solar metal abundance imply a major downward 
revision in Z with respect to the canonical value (i.e., 
from Z Q « 0.02 down to Z Q = 0.0126). If this propa- 
gates to the metallicity scale used for metal-poor stars, 
a major downward revison in the [Fe/H] and Z val- 
ues for Galactic globular cluster stars may be in store 
as well. According to the above discussion, this down- 
ward revision would further complicate the explanation 
of second-parameter pairs in terms of age. The same 
applies, as already stated, if the mass loss results from 
Origlia et al. (2002) are used, as opposed to the mass 
loss recipes given in Table 1 (see Sect. 5 above for a 
critical discussion). 

To close, we note that Grundahl (2003) has recently 
investigated the ages of NGC 288 and NGC 362 on 
the basis of Strbmgrcn photometry, finding that the 
two clusters differ in age by less than 1 Gyr. If so, 
Figure 10 clearly shows that age cannot be the (sole) 
second parameter for this pair. 

7.3.2 The Pair M13-M3 

This pair was first seriously and quantitatively con- 
sidered in the context of the second-parameter phe- 
nomenon by Catelan & de Freitas Pacheco (1995), who 
suggested that the age difference required to account for 
the difference in HB type between M3 (uniformly pop- 
ulated HB) and M13 (blue HB with a long blue tail) 
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Fig. 10 The relative age that is required to explain the 
difference in HB type between NGC 288 and NGC 362 en- 
tirely in terms of age (lines) is plotted as a function of 
the NGC 288 absolute age for two different metallicities: 
Z = 0.001 (upper panel) and Z = 0.002 (lower panel). For 
each panel, each line corresponds to a different mass loss 
formula, as indicated in the insets. The shaded area corre- 
sponds to the turnoff age difference for the pair, from Bel- 
lazzini et al. (2001). (This plot represents an update with 
respect to the similar one presented in Catelan et al. 2001a.) 

exceeded the constraints from turnoff stars. A similar 
conclusion was later reached by Ferraro et al. (1997). 
More recently, Rey et al. (2001) have reanalyzed the 
problem, incorporating the Reimers (1975a, b) mass loss 
"law" into their analysis, and finding an age difference 
of 1.7 ± 0.7 Gyr between the two clusters — which, ac- 
cording to them, "can produce the difference in HB 
morphology between the clusters." Here we revisit the 
problem, investigating the impact of the several differ- 
ent mass loss formulae for red giant stars summarized 
in Sect. 5. 

First wc proceed to compute synthetic HB models 
using the global sample of M3 HB stars discussed in 
Catelan et al. (2001b) and Catelan (2004a), and as- 
suming the canonical mctallicity for the pair — namely, 
Z = 0.001. For M13, wc have retrieved the HST pho- 



Fig. 11 Same as in the previous plot, but for the pair MIS- 
MS. The upper panel shows the comparison between the 
theoretical results and the Salaris & Weiss (2002) turnoff age 
difference, whereas the lower plot compares the theoretical 
values with the Rey et al. (2001) result based on the turnoff 
points. 

tometry from Piotto et al. (2002) and performed num- 
ber counts along the HB of the cluster. When we at- 
tempted to reproduce the observed M13 HB morphol- 
ogy in terms of canonical synthetic HBs, wc quickly 
arrived at the conclusion that a unimodal distribution 
was inadequate, since it was unable to explain the large 
number of very low-mass stars at the extreme hot end of 
the M13 HB. Therefore, we incorporated a second, low- 
mass mode to our simulations, thus obtaining a much 
better agreement with the overall HB morphology of 
the cluster. R should be noted that the EHB of the 
cluster is extremely populous, containing about 30% of 
all HB stars in M13. We then computed the average to- 
tal mass loss required to explain the derived masses of 
the HB stars for each mass loss formula in Table 1 and 
over a range in ages, which then allowed us to compute 
the age difference that was needed, with respect to M3, 
to account for the observed difference in HB type. 
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Fig. 12 Same as in the previous plot, but now ignoring 
M13's EHB component in the calculations. 

The result is shown in Figure 11. In the upper plot, 
we compare the theoretical results with the empirical 
age difference for the pair, as recently derived by Salaris 
& Weiss (2002). (Unless otherwise stated, in this sec- 
tion we use the Salaris & Weiss results obtained in the 
Zinn & West 1984 scale.) In the lower plot, the same 
theoretical results are compared with the age difference 
derived by Rey et al. (2001). We have also compared 
the model results with the VandenBerg (2000) age dif- 
ference values, but these are intermediate between the 
other two studies so that we omit further discussion 
of this case in what follows. Note, in addition, that, 
within the errors, Johnson & Bolte (1998) and Rosen- 
berg et al. (1999) find M13 and M3 to be essentially 
coeval, whereas Stetson (1998) claims that M13 may 
even have the younger turnoff of the two (see below). 
Grundahl (1999) provides an impressive comparison be- 
tween the two clusters in the Stromgren u, (u— y)o plane 
(his Fig. 2, left panel), where one can clearly see that 
there is not much room for a significant age difference 
between the two clusters — Grundahl himself finding an 
age difference of only 0.7 ± 0.2 Gyr between Ml 3 and 
M3, the former being older. 



As can clearly be seen, the required age difference ap- 
pears too large in comparison with the Salaris & Weiss 
(2002) age difference measurements, irrespective of the 
mass loss formula used. The larger age difference mea- 
sured by Rey et al. (2001), on the other hand, may be 
compatible with the turnoff age difference, particularly 
for lower M3 ages and if the "Goldberg" formula pro- 
vides a less reliable description of mass loss rates on the 
RGB than do the others. 

Some authors have suggested, on the other hand, 
that the origin of EHB stars may not be directly linked 
to the age of a globular cluster, being more likely in- 
stead to be due to other physical processes, including 
helium enrichment, binarity, and any other processes 
that may trigger enhanced mass loss on the RGB (e.g., 
Green et al. 2001; Brown et al. 2001). If so, it follows 
that at least 30% of the HB stars in M13 do not owe 
their present-day color to age, a different second param- 
eter being required. In fact, Rey et al. (2001) favor this 
option, although this is in conflict with the scenario of 
Park & Lee (1997) and Lee et al. (2002) for the origin of 
the "UV upturn phenomenon" and for the photomet- 
ric evolution of galaxies, according to which even the 
hottest HB stars owe their existence to age (high ages 
naturally being implied for giant elliptical galaxies in 
this case). Note that this age scenario forms the basis 
upon which the recently measured GALEX ultraviolet 
spectra of early-type galaxies and extragalactic globular 
clusters is currently being interpreted (e.g., Lee et al. 
2005a; Rey et al. 2005). On the other hand, in our own 
galaxy EHB stars are now known to be present even in 
open clusters (Kaluzny & Udalski 1992; Liebert, Saffer, 
& Green 1994; Green et al. 1997), again reinforcing the 
impression that old ages cannot be solely responsible for 
the origin of EHB stars — and hence suggesting that age 
evolution alone cannot fully account for the UV upturn 
phenomenon and for the evolution of the photometric 
properties of galaxies as a function of redshift. 

To check whether the pair M13/M3 might be ex- 
plained by assuming a completely different formation 
channel for the M13 EHB stars, we have repeated the 
above exercise but now removing all M13 EHB stars 
from the sample. Since there are so many EHB stars in 
the cluster, this clearly leads to a significantly higher 
mean mass for the Ml 3 HB stars, and therefore to a 
smaller expected age difference between this cluster and 
M3. Indeed, Figure 12 confirms that, if one admits 
that the EHB stars owe their origin to a different phys- 
ical process, one is able to account for the different HB 
types of M13 and M3 entirely in terms of the reported 
turnoff age differences, irrespective of the adopted mass 
loss recipe, provided (in the Salaris & Weiss 2002 case) 
M3 is younger than 12 Gyr. More stringent constraints 
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on the M3 absolute age would derive in the case of us- 
ing the Origlia et al. (2002) results for the (lack of a) 
dependence of RGB mass loss rate on L, g, or R. 

Finally, it is important to note that, according to 
Stetson's (1998) report on ultra-precise photometry for 
M13 and M3 obtained with the CFHT, the intrinsic 
position of the turnoff point in Ml 3 is actually bluer 
than in MS, which is of course completely unexpected 
if the former is indeed older than the latter. The au- 
thor strongly argues that there are very few (if any) 
sources of systematic errors that could have led him to 
underestimate the turnoff color for M13 compared to 
M3 in his study. One might suspect reddening uncer- 
tainties to be the culprit, but it should be noted that 
both clusters have very low reddening: the Feb. 2003 
edition of the Harris (1996) catalog lists reddening val- 
ues E(B-V) = 0.02 mag and 0.01 mag for M13 and 
M3, respectively (the same values as used by Stetson). 
From Stetson's study, one finds that, for both clusters 
to have the same turnoff color (and hence the same 
ages), the relative reddenings of the two clusters, com- 
pared with the Harris catalog value, would have to be 
incorrect by AE(B — I) = 0.04 mag, which amounts 
to AE(B-V) = 0.015 mag according to the Rieke & 
Lebofsky (1985) standard extinction law — in the sense 
that the M13 reddening must have been overestimated, 
and/or the M3 reddening underestimated. While this 
may seem like a small change, it is worth noting that, 
on the basis of the Schlegel, Finkbciner, & Davis (1998) 
dust maps, one finds E(B-V) = 0.017 mag for M13, 
and E(B-V) = 0.013 mag for M3: while the shifts 
are correct in sign, they are clearly insufficient in size 
to bring M13 to even the same turnoff age as M3. An 
even larger change would be needed, of course, to make 
M13 significantly older than M3. Stetson argues, in 
fact, that a more likely explanation for the differences 
in CMD positions and shapes between the two clusters 
would be provided by a difference in helium abundance, 
as has indeed been suggested by other authors as well 
(e.g., Caloi & D'Antona 2005; Cho et al. 2005). Even 
stronger constraints on the relative reddening values be- 
tween the two clusters would certainly prove helpful in 
clarifying the situation. 

7.3.3 Red HB Globular Clusters in the Extreme Outer 
Halo 

It is worth revisiting the case of the outer halo glob- 
ular clusters with predominantly red HB types which 
have had their ages measured with HST, such as Palo- 
mar 3, Palomar 4, and Eridanus. According to Stet- 
son et al. (1999) and VandenBerg (2000), these clusters 
are slightly younger than the bulk of Galactic globular 



clusters with similar metallicity. Is the detected age dif- 
ference between these clusters and inner halo clusters, 
such as M3 and M5, 19 consistent with the observed dif- 
ference in HB morphology between them? 

The Pair M3-Pal 3 

In terms of metallicity, Pal 3 provides a good match to 
M3. Accordingly, Catelan et al. (2001b) performed a 
study of the age difference between Pal 3 and M3 that 
is required to account for their relative HB types, and 
found that the age difference based on analyses of the 
turnoff points could easily account for the difference 
in mean HB color between the two clusters. We have 
repeated their exercise, finding a slightly larger age dif- 
ference being needed to account for their different HB 
types, but basically confirming their results. On the 
other hand, Catelan et al. call attention to the fact 
that the mass dispersion along the HB of Pal 3 appears 
entirely consistent with zero, thus being significantly 
different from M3's. If differences in mass loss among 
individual red giants is responsible for the presence of 
mass dispersion along the HBs of globular clusters, then 
the mass loss process clearly operated in a different way 
in Pal 3 than it did in M3. Accordingly, Catelan et al. 
(2001a) suggest that while age may be the "global" sec- 
ond parameter driving the mean HB color for at least 
some globular clusters, for many globular clusters envi- 
ronmental effects might be a "local" second parameter 
responsible for generating a dispersion in color around 
this mean value. 

The "Pair" M5-Pal 4/Eridanus 

In terms of chemical composition, both Pal 4 and Eri- 
danus appear to have metallicities similar to M5's. In 
addition, since Pal 4 and Eridanus appear to have sim- 
ilar chemical composition and CMD morphology but 
coarsely populated CMDs, it seems reasonable to com- 
bine the data for the two and perform a single analysis 
of the "pair" comprised by M5, on the one hand, and 
Pal 4/Eridanus, on the other (Catelan 2000). 

Again, we repeat the analysis carried out by Catelan 
(2000), but now comparing the theoretical results with 
the age differences derived on the basis of the turnoff 
points by Stetson et al. (1999), VandenBerg (2000), 
and Salaris & Weiss (2002). The results are shown in 
Figure 13, where the upper panel shows the compar- 
ison between the HB morphology-based analysis and 



19 In fact, proper motion studies (Scholz et al. 1996; Cudworth 
1997) indicate that M5 is an outer halo globular which just hap- 
pens to be close to its perigalacticon at this point in time, the 
cluster actually spending much of its life at galactocentric dis- 
tances larger than 50 kpc. 
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Fig. 13 Same as in the previous plots, but for the case of 
M5 and the outer-halo globulars Pal 4 and Eridanus. The 
latter are studied as a single object, as explained in detail by 
Catelan (2000). The upper panel compares the model pre- 
dictions with the turnoff-based age difference from Stetson 
ct al. (1999) and VandenBerg (2000), whereas the bottom 
panel shows a similar comparison against the turnoff ages 
by Salaris & Weiss (2002). 

the turnoff age difference from Stetson ct al. and Van- 
denBerg, whereas the lower panel shows a comparison 
between the same results for the HB stars in these clus- 
ters and the turnoff ages from Salaris & Weiss. 

These plots (which again give slightly higher HB 
type-based age differences between the clusters than 
originally reported by Catelan 2000) reveal that the 
turnoff age difference measured by Stetson et al. (1999) 
and VandenBerg (2000) is insufficient to account for the 
difference in HB types between M5 and Pal 4/Eridanus, 
irrespective of the mass loss formula used. On the other 
hand, the turnoff age difference reported by Salaris & 
Weiss is clearly more consistent with the hypothesis 
that the age difference between M5 and Pal 4/Eridanus 
is sufficient to account for their relative HB types. It 
is clearly very important to establish what the actual 
turnoff age difference between M5 and Pal 4/Eridanus 



is before we are in a position to decide whether age can 
be the ("global") second parameter for this pair. 

7.4 The Second-Parameter Effect at High Mctallicities 

While the second-parameter problem is traditionally 
thought to affect mostly intermediate-metallicity glob- 
ular clusters, the fairly recent discovery of large (and 
peculiarly bright) RR Lyrae populations (Layden et al. 
1999; Pritzl et al. 2000) and prominent blue HB tails 
(Piotto et al. 1997; Rich et al. 1997) in the moder- 
ately metal-rich Galactic globular clusters NGC 6388 
and NGC 6441 has brought the phenomenon to the 
realm of [Fe/H] ~ —0.5 globular clusters as well. 

Several hypotheses have been discussed in the lit- 
erature to explain the observed HB morphology and 
peculiar RR Lyrae periods in these clusters. These in- 
clude tidal collisions (Rich et al. 1997), helium mix- 
ing on the RGB, a primordially increased helium abun- 
dance, and increased core masses at the RGB tip as a 
result of internal rotation (Sweigart & Catelan 1998). 
Other explanations include a large spread in metallici- 
ties (Piotto et al. 1997; Sweigart 2002), a selective metal 
depletion scenario (Sweigart 1999), and a range in in- 
ternal helium abundances (D'Antona & Caloi 2004). 
In addition, and also as an attempt to explain their 
peculiar HB morphologies, Ree et al. (2002) have sug- 
gested that NGC 6388 and NGC 6441 might be similar 
to u> Ccn in nature, with a (small) internal metallic- 
ity spread and a (fairly large) internal range in ages; 
in their scenario, the blue HB and bright RR Lyrae 
components of NGC 6388 and NGC 6441 would be as- 
cribed to a combination of lower metallicity and old 
ages, RR Lyrae stars being somewhat more metal-poor 
stars evolved away from a position on the blue ZAHB. 

Unfortunately, most — if not all — of these scenarios 
face strict observational and theoretical constraints. 
Rich et al. (1997) show that tidal collisions cannot pro- 
duce the observed HB morphology — and we add that it 
cannot lead to peculiarly bright RR Lyrae stars, either. 
A primordially increased helium abundance does not 
appear consistent with the position of the RGB bump 
in these clusters (Raimondo et al. 2002). A large spread 
in metallicities, as also pointed out by Raimondo et al., 
is not supported by observations of the cluster CMDs 
in the RGB region, since the latter do not show the 
large spread in colors that would be expected in this 
case. Moreover, Clementini et al. (2005) have recently 
found, based on VLT spectra, only a small (though sig- 
nificant) spread in [Fe/H] among the RR Lyrae stars in 
NGC 6441 — although this result could not be confirmed 
by Gratton ct al. (2007) among the cluster's red giants. 
Evolution away from a position on the blue ZAHB does 
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not produce enough bright RR Lyrae variables to ex- 
plain the observed RR Lyrae period distributions, and 
neither is the sloping nature of the HB quantitatively 
reproduced in the Ree et al. (2002) scenario (Pritzl et 
al. 2002b). 

Most puzzling of all is the evidence that the blue 
HB stars in these clusters cooler than the Grundahl et 
al. (1999) jump do not appear to have peculiarly low 
gravities (Moehler, Swcigart, & Catclan 1999a), which 
effectively rules out any scenario that requires the blue 
HB and RR Lyrae components to be brighter than in 
canonical models. However, a brighter HB at both the 
blue HB and the RR Lyrae levels seem to be required 
by i) The glaring evidence that the tip of the blue HB 
is ~ 0.5 mag brighter in V than the red HB compo- 
nent; ii) The very long periods of the RR Lyrae stars 
in both clusters. Given the apparently irreconcilable 
evidence, we suggest that a reassessment of the spec- 
troscopic gravities for a larger sample of blue HB stars 
in both clusters would prove well worth the effort. 20 

Assuming that the spectroscopic gravities can indeed 
be reconciled with the requirement that the blue HB 
component be unusually bright, we discuss the pos- 
sibility that a internal spread in helium abundances 
(D'Antona & Caloi 2004) may help account for the 
HB morphology and RR Lyrae pulsation properties in 
these clusters (see Sweigart & Catelan 1998 for numer- 
ical simulations in the case of helium mixing) . 

Consider the case of NGC 2808, and assume that, 
as suggested by D'Antona & Caloi (2004), there is an 
internal helium abundance range among the stars in 
the cluster. In their scenario, NGC 2808's red HB and 
RR Lyrae components have a "normal" helium abun- 
dance, blue HB stars at the "horizontal" level have a 
mild helium enhancement, and hotter blue HB stars 
have a much higher (initial) helium content, up to about 
35%. Lee et al. (2005b) suggest a similar scenario for 
the cluster, only invoking an even higher helium en- 
hancement for the blue HB stars. Now imagine how 
NGC 2808 might have looked like in the not too remote 
past (say, a couple of Gyr ago): its present-day helium- 
enhanced blue HB stars would necessarily be redder as 
a consequence of the younger age, and therefore helium- 
enhanced stars would be present in relatively high num- 
bers inside the RR Lyrae instability strip — thus leading 
to overluminous, long-period RR Lyrae stars. The stars 
at the tip of the blue HB would have an even higher he- 
lium abundance, leading to a marked sloping HB when 
compared with the red HB and RR Lyrae components 



20 After this paper had been completed, Moehler & Sweigart 
(2006) revisited the problem and concluded that those previous 
results were indeed incorrect, most likely due to problems with 
the background subtraction from the spectra (blends). 



in the same cluster. Clearly, the NGC 2808 envisaged 
by D'Antona & Caloi and Lee et al. would have looked 
a lot like the present-day NGC 6388 and NGC 6441 in 
the past! Accordingly, a possible explanation for the ob- 
served peculiarities in these clusters could involve both 
an enhanced helium abundance among a fraction of the 
cluster stars (note that the blue HB plus RR Lyrae com- 
ponents in both clusters constitute relatively small frac- 
tions of the overall HB populations, so that the RGB 
bump constraint might not be violated in this case) and 
a younger age with respect to the bulk of the Galactic 
globular clusters — the latter not necessarily being re- 
quired, depending on mctallicity (i.e., first parameter) 
and mass loss effects: indeed, the deep HST photom- 
etry presented by Catelan et al. (2006) for NGC 6388 
reveals that the cluster is comparable in age to 47 Tuc. 

A potential problem with this scenario is the fact 
that, for every single helium abundance value, one 
should expect, by analogy with what is observed in 
"single-population globular clusters," a spread in mass 
loss as well. Therefore, at any given color, a spread in 
helium abundance should also be present, and it is un- 
clear whether the implied HB luminosity distribution 
would match in detail the observed one, which appears 
fairly tight. This problem is especially evident in the 
recent simulated CMD by Busso et al. (2007; see their 
Fig. 7), but is not as apparent in the simulated CMDs 
of Caloi & D'Antona (2007) — which again may be due 
to the lack of a mass spread for the populations with 
enhanced Y in their simulations. Note, in this sense, 
that the simulations presented by Sweigart & Catelan 
(1998) in their helium mixing and rotation scenarios 
also ignore the effect of a spread in masses for each 
individual Y. In any case, Piotto (2008) has recently 
reported that NGC 6388 does show a double subgiant 
branch, thus reinforcing the evidence for the presence 
of more than one stellar population in the cluster. 

We conclude that stringent tests of such a scenario 
could be provided by deep HST photometry (so that 
the turnoff ages, and possible splits indicative of mul- 
tiple stellar populations, can be reliably established), 
spectroscopic gravities for a large sample of moderately 
cool blue HB stars in both clusters, and HB simulations 
in which both a spread in Y and in RGB mass loss are 
simultaneously taken into account. 

8 On the RR Lyrae Luminosity-Metallicity 
Relation 

Much has been written over the past several years in 
regard to the luminosities of HB stars at the RR Lyrae 
level, particularly in the V band, and its dependence on 
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metallicity; recent reviews of the subject include papers 
by Chaboyer (1999), Cacciari (1999, 2003), Cacciari & 
Clementini (2003), and Storm (2006). Bono (2003) and 
Cassisi (2005) have provided recent reviews in which 
theoretical uncertainties affecting the predicted prop- 
erties of HB stars, including their luminosities, have 
been discussed in considerable detail. Several different 
theoretical results have recently been compared by Cac- 
ciari & Clementini (2003; see their Fig. 1) and Gallart, 
Zoccali, & Aparicio (2005; see their Fig. 9). Accord- 
ingly, in what follows we will content ourselves with 
providing an updated estimate of the HB luminosity- 
metallicity relation as based on the trigonometric par- 
allax of RR Lyrac itself, only briefly mentioning recent 
progress in understanding some discrepancies that have 
prevented the establishment of a universally accepted 
My(RRL) - [Fe/H] relation. 

It should be noted that this crucial relation has tra- 
ditionally provided the very basis for the Population II 
distance scale, thereby constituting one of the most im- 
portant techniques used to help nail down the first step 
in the cosmological distance ladder — namely, the dis- 
tance to the LMC (e.g., Benedict ct al. 2002; Kovacs 
2003; Storm 2006). Moreover, its slope and zero point 
have long been recognized as crucial ingredients in the 
determination of the absolute ages of globular clusters 
and their variation with metallicity (e.g., Sandage & 
Cacciari 1990; Walker 1992c). On the other hand, it 
should also be noted that this relation can only be used 
in a very approximate way to estimate the average ab- 
solute magnitudes of RR Lyrae stars of a given metal- 
licity. 

In order to properly evaluate the absolute magni- 
tudes of individual RR Lyrae stars, more precise tech- 
niques, frequently involving a period-luminosity rela- 
tion, are required. Unlike the case of classical Cepheids, 
however, good period-luminosity relations are not avail- 
able for RR Lyrae stars in the visual bandpasses, for 
reasons that have been discussed in detail by Catelan, 
Pritzl, & Smith (2004). On the other hand, it has long 
been known that good RR Lyrac period-luminosity re- 
lationships are present in the near-infrared (Longmore, 
Fernley, & Jameson 1986; Longmore et al. 1990). Em- 
pirical results have recently been critically discussed by 
Sollima, Cacciari, & Valenti (2006) and Feast et al. 
(2008), and Sollima et al. (2008) have recently provided 
a detailed analysis of the J, H , K light curves of RR 
Lyrae itself. Theoretical calibrations using the near- 
infrared bandpasses J, H, K have been provided by 
Cassisi et al. (2004), Catelan et al., and Del Principe et 
al. (2006), among others. As pointed out by Soszynski 
et al. (2003) and Catelan et al., a period- luminosity re- 
lation is also present in /, although in this case one 



finds somewhat more scatter and a stronger metallicity 
dependence than in the near infrared. Finally, Cortes 
& Catelan (2008) and Caceres & Catelan (2008) have 
recently shown that very precise period-luminosity and 
period-color relations may also be defined for RR Lyrae 
stars in the Stromgren and SDSS filter systems, respec- 
tively. 

8.1 The Variable Star RR Lyr and the HB 
Luminosity-Metallicity Relation 

The star RR Lyr is the closest of its class, and accord- 
ingly has proven of great interest for trigonometric par- 
allax studies. Its variability was noted by Williamina 
P. Fleming at the Harvard College Observatory prior to 
July 1889, but the discovery was not announced until 
a few years later by Pickering (1901). 

Benedict ct al. (2002) obtained, using the Hubble 
Space Telescope, a much more accurate (and signifi- 
cantly smaller) value for the absolute parallax of RR 
Lyrae than had previously been provided by Hipparcos, 
namely 7r a b s — 3.82 ± 0.20 mas (compared to ir^s = 
4.83 ± 0.59 mas; Perryman et al. 1997). Recently, 
van Leeuwen (2007) revised the trigonometric paral- 
laxes provided by Hipparcos, and arrived at a value 
^abs = 3.46 ± 0.64 mas for RR Lyr. A weighted aver- 
age of ground-based studies (van Altena, Lee, & Hof- 
fleit 1995) indicated a parallax 7r a b s = 3.0 ± 1.9 mas for 
RR Lyrae (see Fig. 6 in Benedict et al. 2002). Taking a 
weighted average of these results, we obtain a final value 
of 7r a b s = 3.78 ± 0.19 mas for RR Lyr. This implies a 
revised distance modulus of (to— M) = 7.11 ±0.11 mag 
for the star. 

Benedict et al. (2002) argue in favor of a rela- 
tively low extinction towards RR Lyr, namely Ay ~ 
0.07 ± 0.03 mag. In recent work, an intensity-weighted 
mean magnitude of (V") = 7.76 mag (Fernley et al. 
1998a) 21 has been adopted for RR Lyr. However, we 
note that this value is based on Hipparcos photometry, 
which may require a non-trivial transformation to the 
standard system. For comparison, Layden (1994) de- 
termines a (V) = 7.66 mag, and Layden ct al. (1996) 
find instead (V) = 7.74 mag. On the other hand, 
Gould & Popowski (1998) argue strongly in favor of 
the Hipparcos-b&sed magnitudes of Fernley et al., only 
proposing an additional, reddening-related correction: 
Vgpqq = Vf98 — 0.2 E(B — V). Using a standard ex- 
tinction law, Ay ~ 3.1 E(B — V), and the reddening 
value obtained by Catelan & Cortes (2008) on the 



21 This value is provided in Table 1 of their paper, 
which is only available in electronic format, from 

http: //vizier . u-strasbg. fr/viz-bin/VizieR?-source=J/A+ A/330/515. 
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basis of Stromgren photometry, namely E(B — V) = 
0.015 ±0.020 mag, leads to a final, extinction-corrected 
RR Lyr mean magnitude of (Vb} — 7.71 ± 0.06 mag. 

It is important to note that the intensity-mean mag- 
nitude (and even more so the corresponding magnitude- 
weighted average) does not necessarily correspond to 
the magnitude of the "equivalent static star" (i.e., the 
magnitude the star would have if it were not pulsat- 
ing): an amplitude-dependent correction has to be 
applied. Such a correction has been obtained, both 
for fundamental (RRab) and first-overtone (RRc) vari- 
ables, by Bono, Caputo, & Stellingwerf (1995) on the 
basis of detailed hydrodynamical pulsation models of 
RR Lyrae stars. RR Lyrae has long been known to 
present the Blazhko (1907) effect, and its Peiazhko ap- 
pears to fall in the range between 40 d (Smith et al. 
2003) and 41 d (Szeidl & Kollath 2000), with an ad- 
ditional, longer-term periodicity (P ~ 4 yr) also being 
present (Detre & Szeidl 1973). Inspection of the light 
curves for RR Lyr presented by Smith (1995), Szeidl 
& Kollath, and Smith et al. suggest that the ampli- 
tude in V oscilates in the range between 0.5 mag and 
1.1 mag. According to Table 2 in Bono et al., this is 
precisely the amplitude range over which the intensity- 
weighted mean magnitude provides the most accurate 
description of the magnitude of the equivalente static 
star. Therefore, no amplitude corrections appear to be 
needed in the case of RR Lyrae. Taking, accordingly, 
a value (Vb) = 7.71 ± 0.06 mag and a distance modu- 
lus (to — M)o = 7.11 ± 0.11 mag, one finds an absolute 
magnitude for the star of My = 0.60 ± 0.13 mag. 

In regard to the star's metallicity, Clementini et 
al. (1995) obtain [Fe/H] = -1.39 dex and [a/Fe] = 
0.31 dex, which are quite typical values for Galac- 
tic halo stars. Previous metallicity measurements for 
this star had provided values in the range between 
[Fe/H] = -1.14 and -1.21 dex (see Table 6 in Clemen- 
tini et al. 1995), so that the Clementini et al. result rep- 
resented a significant downward revision. According to 
Bragaglia et al. (2001), the Clementini et al. measure- 
ments are in a scale that closely mimicks the Zinn & 
West (1984) scale. We accordingly adopt a metallicity 
value [Fe/H] = -1.39 ± 0.10 dex for RR Lyr. 

Note that the analysis of Catelan & Cortes (2008) 
suggests that RR Lyr is a somewhat evolved star, with 
an ovcrluminosity of ~ 0.06 ± 0.01 mag in V with re- 
spect to the mean for other RR Lyrae stars with similar 
metallicity (see also Feast et al. 2008). To within the 
errors, this result is not inconsistent with the recent 
near-infrared study by Sollima et al. (2008). There- 
fore, the average luminosity of RR Lyrae variables with 
metallicity similar to that of RR Lyr itself should be 
around 0.66 ± 0.13 mag. 
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What do these results imply, in terms of the tra- 
ditionally employed RR Lyrae luminosity-metallicity 
relation in the V band, usually taken in the form 
My(RRL) = a [Fc/H] + (31 

To answer this question, we shall first adopt a slope 
a = 0.23 ± 0.04 for this relation, as found and/or fa- 
vored in several recent reviews of the subject, including 
Chaboyer (1999), Cacciari (1999, 2003), and Cacciari 
& Clementini (2003). Several recent analyses do pro- 
vide additional support for this result: for instance, 
Clementini et al. (2003) and Gratton et al. (2004a) ob- 
tain a = 0.214 ± 0.047 from analysis of RR Lyrae vari- 
ables in the LMC, whereas Olech et al. (2003) obtained 
a = 0.21 — 0.28, depending on their treatment of pre- 
sumably well-evolved RR Lyrae variables with periods 
around 0.7 d, from analysis of the RRab stars in u) Cen. 

Using the Clementini et al. (1995) metallicity for RR 
Lyr in the Zinn & West (1984) scale, one then finds a 
value of P = 0.98 ± 0.13, implying a final relationship 
of the following form: 

My(RRL) = (0.23±0.04) [Fc/H] zw + (0.98±0.13). (3) 

If one transforms the Clementini et al. metallicity to 
the Carretta & Gratton (1997) scale and then repeats 
the analysis, one finds instead 

My(RRL) = (0.23 ±0.04) [Fc/H] CG + (0.93 ±0.13). (4) 

These relations, while based on a detailed reassess- 
ment of the absolute magnitude and evolutionary sta- 
tus of RR Lyr, turn out to be similar to the rela- 
tion derived in the recent review papers by Chaboyer 
(1999), Cacciari (1999, 2003), and Cacciari & Clemen- 
tini (2003), based on a critical analysis of several cali- 
bration techniques (but ignoring the evolutionary sta- 
tus of the star). This notwithstanding, some meth- 
ods have provided somewhat discrepant slopes and/or 
zero points, and we shall momentarily address two such 
cases. Before doing so, however, we will immediately 
proceed to deriving the all-important distance modulus 
of the LMC that is implied by these relations (see also 
Alves 2004 for a recent review). 

8.2 The Distance Modulus of the LMC 

Eq. (3) implies an absolute magnitude My(RRL) = 
0.64 ± 0.14 mag at [Fe/H] = -1.48 ± 0.07. The lat- 
ter is the mean metallicity derived for LMC RR Lyrae 
variables by Gratton et al. (2004a), in the Zinn & West 
(1984) scale. Using a value (V ) = 19.068 ± 0.102 mag 
from Gratton et al., one then finds an updated true dis- 
tance modulus for the LMC of {m-M)^ MC = 18.42 ± 



34 



M. Catelan 



0.17. If one uses instead the average values for LMC 
RR Lyrae stars independently determined by Borissova 
et al. (2004), namely [Fe/H] = -1.46 ± 0.09 dex and 
(V) = 19.45 ±0.04 mag, with their favored reddening of 
E(B-V) =0.11 mag, one finds (V ) = 19.11±0.04 mag 
for a My = 0.65 ± 0.14 mag, thus implying a true dis- 
tance modulus (m-M)^ MC = 18.46 ± 0.15. Taking a 
weighted average over these two results, we arrive at 
the following distance modulus for the LMC, based on 
our updated analysis of the star RR Lyr: 

(m-M)^ MC = 18.44 ± 0.11. (5) 

8.3 Are We Converging on a My (RRL) - [Fe/H] 
Relation Yet? 

In what follows, we discuss a few discrepant calibrations 
of the HB luminosity-metallicity relation, and describe 
how the problem has recently been solved or what sug- 
gestions may have been advanced to reconcile the dis- 
crepant calibrations. 

8.3.1 The Shallow Slope Obtained from HST 
Photometry of M31 Globular Clusters 

Almost a decade ago, a very shallow slope, a = 0.13 ± 

0. 07, was derived by Fusi Pecci et al. (1996) on the ba- 
sis of HST observations of M31 globular clusters. How- 
ever, the determination of the HB level in their CMDs 
was far from straightforward, since instead of seeing a 
horizontal branch at the RR Lyrae level, whenever a 
blue or intermediate HB component was present their 
CMDs revealed instead surprisingly sloped "horizontal" 
branches, not unlike what one sees when plotting an 

1, (V-T) CMD for Galactic globulars. Such sloping HBs 
are not seen among Galactic globular clusters, and even 
though differences between Galactic and extragalactic 
globulars (related, for instance, to the different chem- 
ical evolution histories of the different galaxies, or to 
the amount of angular momentum available in the dif- 
ferent protogalactic clouds) may indeed exist, theoret- 
ical models of HB stars do not predict similar CMD 
morphologies. 

A possible solution to this puzzle has recently been 
provided by Rich et al. (2005), who found that the Fusi 
Pecci et al. (1996) CMDs were strongly affected by pho- 
tometric blends, not accounted for in their original anal- 
ysis. Investigating the impact of these blends upon the 
morphology of the HB in their CMDs, Rich et al. have 
found that the strange sloping nature of the HB in Fusi 
Pecci et al. is likely due to the presence of photometric 
blends. Accordingly, Rich et al. derived new CMDs 



for an enlarged sample of M31 GCs in which the ef- 
fects of blends were taken into account. As a result, 
they provide a revised slope of a ~ 0.20 ± 0.09, clearly 
in much better agreement with eq. (3). It should be 
noted, however, that careful inspection of the CMDs 
published by these authors still reveal unrealistic HB 
shapes, thus raising the possibility that additional cor- 
rections will be needed before a final relation between 
HB magnitude and metallicity can be derived on the 
basis of observations of M31 globulars. 

8.3.2 The Faint Zero Point of the Method of 
Statistical Parallaxes 

As far as the zero point of eq. (3) is concerned, the 
recent review by Cacciari & Clementini (2003) shows 
that there is reasonable agreement among the several 
different methods that are used to infer it. Impor- 
tantly, the Baade-Wesselink method, which used to fa- 
vor a faint zero point (i.e., fainter than provided by 
the above calibration by at least 0.1 mag; see review 
of earlier work by Fernley et al. 1998b), is now seen to 
be in good agreement with eq. (3) and with a brighter 
zero point for the HB luminosity-metallicity relation. 
Indeed, Kovacs (2003), by applying the same Baade- 
Wesselink algorithm to both Galactic RR Lyrae and 
Cepheid variables, has recently shown that a consistent, 
"long" distance modulus for the LMC obtains, namely 
(m — M)o = 18.55 mag — about 0.05 mag brighter than 
implied by eq. (3). 

This notwithstanding, at least one method remains 
that does keep repeatedly providing a faint zero point 
to the LMC: the statistical parallaxes method. For 
instance, Gould & Popowski (1998) favor a My = 
0.77 mag at [Fe/H] = —1.6 dex from analysis of a 
sample of 147 RR Lyrae stars, which is 0.21 mag 
fainter than implied by eq. 3; or My = 0.80 mag at 
[Fe/H] = —1.7 dex after adding to the sample 716 
non-variable metal-poor stars, which is 0.26 mag fainter 
than eq. (3). While Cacciari & Clementini (2003) hint 
that the presence of disk stars may affect the Gould & 
Popowski results, Dambis & Rastorguev (2001) recently 
provided a new application of the method in which thick 
disk stars were carefully separated from halo stars using 
both kinematic and metallicity criteria, but essentially 
confirming the Gould & Popowski results. Popowski & 
Gould (1998a, b) have provided a very careful and de- 
tailed analysis of the possible systematic uncertainties 
affecting the method, without succeeding in identify- 
ing a likely cause for the difference with respect to sev- 
eral other methods — in fact, they hint that the problem 
lies with the latter (see also Popowski & Gould 1999). 
It is indeed unclear what the solution to this problem 
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will be, but Cacciari & Clementini suggest that the 
Dambis & Rastorguev results may be affected by inho- 
mogeneities in the distribution of their halo stars. In- 
deed, that the distribution of Galactic halo stars is not 
quite uniform has recently been noted, using blue HB 
stars, by Altmann et al. (2005) and Clcwlcy ct al. (2005) 
(but see also Brown et al. 2004a); using RR Lyrae stars, 
by Ivezic et al. (2000), Vivas et al. (2001), Vivas & Zinn 
(2003, 2006), Kinman et al. (2004), Zinn et al. (2004), 
and Keller et al. (2008); and, using several other types 
of tracers, by, for instance, Newberg et al. (2002), Ibata 
et al. (2003), Majewski et al. (2003, 2004), Martmez- 
Delgado et al. (2004), Newberg & Yanny (2005), and 
Rocha-Pinto et al. (2006), among many others. 22 

8.3.3 Differences between HB Stars in Globular 
Clusters and in the Field? 

The similarity between metal-poor field and cluster RR 
Lyrae stars appears reasonably established (Catelan 
1998; De Santis & Cassisi 1999; Carretta et al. 2000), 
but there may be differences in regard to at least the 
more extreme HB stars, both in regard to their luminos- 
ity distribution (Brown et al. 2005 claim that the field 
contains significantly fewer stars at the blue end than 
do such globular clusters as M15 and NGC 6229) and 
their physical origin: Moni Bidin ct al. (2006a, b) and 
Moni Bidin, Catelan, & Altmann (2008) find that the 
EHB stars in NGC 6752 and M80 (NGC 6093) do not 
appear to be associated to binary systems, unlike what 
seems to happen most frequently among field sdB stars 
(Green et al. 2001). Moni Bidin et al. (2008) interpret 
this phenomenon as likely due to a binary fraction-age 
relation (globular cluster stars being on average much 
older than field stars), a hypothesis which has recently 
been supported by Han (2008). In addition, the smaller 
envelope masses of the older cluster red giants make it 
much easier for the single-star channel to produce EHB 
stars, as also noted by Moni Bidin et al. (2008). 

It should also be noted that the abundance anoma- 
lies which are commonly seen in globular cluster red 
giants (e.g., Sneden ct al. 2004; Gratton et al. 2004b; 
and references therein) appear not to be present in field 
red giants (Gratton ct al. 2000; see also Sect. 5.2 in 
Grundahl et al. 1999 for extensive references to work 
prior to the year 2000). Therefore, if these abundance 
anomalies affect somehow the evolution of RGB and HB 



22 The widespread use of HB stars, whether variable or not, in 
many such studies clearly shows their importance as probes of 
the Galaxy's structure, formation, and evolution. Other recent 
examples include Sirko et al. (2004), Thorn ct al. (2005), Brown 
ct al. (2005, 2008), Kaempf, de Boer, & Altmann (2005), and 
Clcwlcy & Kinman (2006) — among others. 



stars (through the "primordial" and/or "deep mixing" 
channel), as has indeed been often suggested in the lit- 
erature (see §7.1 for extensive references), they should 
naturally be expected to have an impact upon the ob- 
served properties of HB stars — including color distribu- 
tion, luminosities, gravities, and pulsation characteris- 
tics. While these differences have not been frequently 
detected, it is worth recalling the cases of NGC 6388 
and NGC 6441, whose RR Lyrae stars are very differ- 
ent from field RR Lyrae stars with similar metallicity 
(Sweigart & Catelan 1998; Layden et al. 1999; Pritzl et 
al. 2000, 2001, 2002b; Clementini et al. 2005; Catelan 
et al. 2006; Matsunaga et al. 2006; Matsunaga 2007; 
see also Layden 1995 and Catelan 2004b for additional 
references to metal-rich globular clusters harboring rel- 
atively long-period RR Lyrae stars). 

8.34 The My(RRL) - [Fc/H] Relation: Linear, 
Quadratic, or Even More Complex? 

Steeper slopes than provided by eq. (3) may result from 
the inclusion of RR Lyrae stars more metal-rich than 
[Fc/H] ~ —1, as originally suggested by Castellani, Chi- 
effi, & Pulone (1991) from theoretical computations. In 
fact, for field stars in our galaxy, a quadratic relation 
between My (RRL) and [Fe/H] is likely to be superior to 
a linear approximation when such metal-rich variables 
are included, in the sense that the metal-rich variables 
are fainter than would be expected using a linear fit to 
the more metal-poor data — a result which is strongly 
supported by the theoretical models (e.g., Castellani et 
al. 1991; Caputo et al. 2000; Bono ct al. 2003; Catelan 
et al. 2004; and references therein). 

On the other hand, it should also be noted that the 
helium abundance is also expected to increase with Z , 
by amounts which may differ in different environments 
as a consequence of different chemical enrichment laws 
(as early noted by van den Bergh 1967). Accordingly, 
it is extremely important to realize that there is no 
strong physical basis for a universal Y — [Fc/H] relation. 
As an example, Catelan & de Freitas Pachcco (1996) 
analyzed three different chemical enrichment scenarios 
which, though providing essentially the same helium 
abundance at low metallicities ([Fe/H] < —1), led to 
differences in Y by ~ 0.035 at [Fe/H] « -0.5, and by 
up to 0.1 at solar metallicity (see also Catelan 2008b). 
Indeed, it would perhaps be surprising if galactic disks 
and spheroids presented precisely the same AY — AZ 
"law," given the evidence that they present different a- 
element enrichment patterns — variations in the latter 
being expected to be accompanied by variations in the 
AY — AZ relation as well (see Figs. 1 and 2 in Catelan 
2008b). 
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This can impact HB luminosities in a quite dramatic 
way. Since the HB luminosity depends on the helium 
abundance as 



(Catelan 1996), 23 the quoted differences in helium 
abundance may lead to changes in the HB luminos- 
ity, from one chemical enrichment scenario to the next, 
by as much as 0.15 mag at [Fe/H] = —0.5, and by 
0.45 mag (!) at [Fe/H] = 0. The helium enrichment 
law clearly plays a crucial role in defining the slope of 
the HB luminosity-mctallicity relation at high metal- 
licities, and this should be properly taken into account 
when studying different stellar populations. 24 

In the same token, we believe it is a risky proce- 
dure to derive relative globular clusters ages, using the 
HB luminosity level as a standard candle, for globular 
clusters of different metallicities belonging to different 
populations, as has been done by Ortolani et al. (1995) 
in the case of 47 Tuc and NGC 6528, since they may 
have undergone difference helium enrichment laws and 
therefore have significantly different HB luminosities. 

Another important implication of this result is that 
the production of hot HB sources at high metallicities, 
which are believed to be the main sources of the "UV 
upturn" (or UVX) light coming from elliptical galax- 
ies and the bulges of spirals (see O'Connell 1999 for 
a review and extensive references), is importantly af- 
fected by their precise helium enrichment laws (e.g., Yi 
et al. 1997). In the Galactic bulge, evidence of the pres- 
ence of hot HB stars has been provided and discussed 
by Bertelli et al. (1996), Terndrup et al. (1999, 2004), 
Peterson et al. (2001), Busso et al. (2005), and Busso 
& Moehler (2008). It should also be noted that the 
precise RGB mass loss recipe adopted may also signif- 
icantly impact the predicted production of UV sources 
(see Sect. 5). 



23 This is a little higher than what obtains from ZAHB models 
(see, e.g., Table 1 in Swcigart, Rcnzini, & Tornambc 1987), due 
to the fact that high-y evolutionary tracks present much more 
important luminosity evolution than do models for lower Y (see, 
e.g., Sweigart & Gross 1976; Swcigart 1987). 

24 We note, in passing, that possible differences in the helium en- 
richment law should also be taken into account when studying 
classical Cepheids in external galaxies, fn particular, Fiorentino 
et al. (2002) show that the Cepheid distance scale also presents 
an important, intrinsic dependence on Y. Note that the usual 
approach, in Cepheid-based studies of the cxtragalactic distance 
scale, is to take into account statistical uncertainties in a "univer- 
sal" AY — AZ enrichment law, whereas our point is that such a 
relation need not be universal. Fiorentino et al. already point out 
that the galaxy M101 may require different enrichment laws in its 
inner and outer parts in order to account for observed differences 
in their dependence on metallicity. 



Therefore, a calibration such as the one provided 
by eq. (3) may not be straightforwardly applicable to 
just any galaxy, especially when more metal-rich com- 
ponents are considered. Also, from Rittcr's (1879) re- 
lation, one finds that the pulsation properties of RR 
Lyrae stars are directly inherited from, or rather re- 
flected upon, their luminosities. However, Figure 5 
clearly shows that the RR Lyrae that belong to our 
galaxy present a rather peculiar behavior, at least 
compared to those belonging to its neighboring galax- 
ies. But if this is indeed so, what guarantees do we 
have that iWy(RRL) — [Fc/H] relations based upon 
Galactic halo RR Lyrae stars are universally applica- 
ble to other galaxies? Arguably, one should make sure 
that, if not the Oosterhoff dichotomy itself, at least 
the Oosterhoff- Arp-Preston-Sandage period-metallicity 
progression (Oosterhoff 1939, 1944; Arp 1955; Preston 
1959; Sandage 1981, 2004; and references therein) is 
verified in the sample under scrutiny (Catelan 2004b) — 
which is only possible by using time-series observations. 

Even for bona-fidc Galactic RR Lyrae stars, one 
should recall that the My(RRL) — [Fe/H] relation can 
be strongly affected by evolutionary effects (e.g., Lee 
1990; Demarquc et al. 2000); these are more properly 
avoided when using the RR Lyrae period-luminosity re- 
lation to infer distances. In any case, one must keep 
in mind that there appear to be "pathological outliers" 
for which the available calibrations of any of these re- 
lations may not be straightforwardly applied, due to 
the fact that they may not be properly described by 
canonical models (Sweigart & Catelan 1998; Pritzl et 
al. 2001, 2002b; Raimondo et al. 2002). Also, if the 
helium abundance is enhanced in at least a fraction 
of the stars in some globular clusters, whether due to 
primordial or to evolutionary effects, models assuming 
Y ~ 0.23 — 0.25 will not provide a correct distance 
to the RR Lyrae stars in the systems which partook a 
primordial enhancement, or in which some RGB stars 
may have undergone helium enrichment and ended up 
as over luminous RR Lyrae. For critical discussions of 
several different channels for the production of high-Y 
stars in globular clusters, the reader is referred to the 
extensive list of references provided in §7.1. 

Even if present-day Galactic globular clusters have 
only their blue HB components affected by a high Y, 
a few Gyr ago these blue HB components were actu- 
ally (overluminous) RR Lyrae stars. Therefore, even 
if not in our present-day galaxy, there is no guarantee 
that in other (somewhat younger) galaxies and their 
respective globular cluster systems we might not run 
into overluminous RR Lyrae that are their analogues 
of the present-day, helium-enhanced, Galactic blue HB 
stars. The effect of a change in Y upon the RR Lyrae 
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period-luminosity relation has been recently discussed 
by Catelan et al. (2004). Again, over/under luminosity 
should be reflected upon the pulsation periods of RR 
Lyrae stars, so that proper monitoring of their pulsa- 
tion cycles remains the only way to safely apply local 
(or standard theoretical) HB luminosity calibrations to 
external systems. For nonvariable HB stars, unfortu- 
nately such a consistency diagnostic is not available. 

For the sake of argument, consider the Yoon & Lee 
(2002) scenario, whereby the Galactic OoII component 
is not a genuine component of our galaxy, having in- 
stead been accreted from an external galaxy at some 
point in the past. 25 In this scenario, had such an ac- 
cretion event not taken place, our Galactic globular 
cluster system would be characterized by the Ool com- 
ponent shown in Figure 5 (left panel) plus the OoIII 
component shown in the same panel. In this case, one 
would find for our Galactic globulars an opposite trend 
to what is seen in the Oostcrhoff-Arp-Preston-Sandage 
progression, with the more metal-rich globulars pre- 
senting much longer periods — and therefore presumably 
having brighter HBs — than the Ool clusters. There is 
no guarantee that in external galaxies similar to our 
own but which have not (again in the Yoon & Lee sce- 
nario) incorporated an "OoII protogalactic fragment" 
in the course of its history, such a seemingly far-fetched 
behavior is not precisely what happened in practice. 

8.4 A Caveat: The Evolutionary Status of the Star 
RR Lyr 

One possible caveat with employing the star RR Lyr 
to constrain the average HB luminosity-metallicity re- 
lation at the instability strip level is the fact that we do 
not know a priori the evolutionary status of the star. 
Therefore, it could be significantly brighter (or fainter) 
than the average field RR Lyrae at the same metallic- 
ity, in which case our results could be systematically off 
by 0.1 mag or more, given the intrinsic spread in RR 
Lyrae luminosities (e.g., Sandage 1990). Indeed, Cate- 
lan & Cortes (2008) have recently argued, based on 
Stromgren photometry and theoretical modelling, that 
RR Lyr may be brighter than other field RR Lyrae stars 
of similar metallicity by 0.064±0.013 mag in V, whereas 
Feast et al. (2008), comparing distances derived using 
type II Cepheids and RR Lyrae stars, similarly advo- 
cate that RR Lyr is over luminous by 0.16±0.12 mag 
in V. Therefore, it may not be entirely safe to assume, 
as frequently done, that RR Lyr is truly representative 

25 In fact, van den Bcrgh (1993) and Lcc & Carney (1999b), 
among others, suggest precisely the opposite, namely, that it is 
the Ool component of the Galaxy that is of external origin, the 
OoII clusters having formed very early in the Galaxy's infancy. 
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Fig. 14 Observed minus computed (O-C) phase of maxi- 
mum light (in days) as a function of heliocentric Julian day 
for the variable star RR Lyrae. 

of other RR Lyrae-type stars of similar metallicity. A 
definitive solution to this problem will probably have to 
wait for new, accurate parallax determinations for large 
numbers of RR Lyrae stars, such as will be afforded by 
the GAIA and SIM missions. 

In the meantime, however, it could also be useful to 
look further into another diagnostic that is available to 
us at this point in time: the period change rate of RR 
Lyr. More specifically, a large, positive period change 
rate for the star could raise a yellow flag indicating the 
possibility of its being in a fast, advanced evolution- 
ary stage, and therefore likely brighter than most other 
stars of similar metallicity The presence of a near-zero 
or negative period change rate would be more in line 
with the star being found around the slowest part of 
the HB evolution. Of course, the caveat should be kept 
in mind that it is very often the case that RR Lyrae 
stars are seen to exhibit erratic period changes which 
are in some cases even orders of magnitude different 
from what is predicted by stellar evolution theory (see 
Smith 1995 for a review and discussions), so that the 
present argument should be taken with due caution. 

Automatically computed and updated values of ob- 
served minus computed (O-C) times of maxima for RR 
Lyr can be found in the GEOS RR Lyr Database web 
page 26 (see Le Borgne, Klotz, & Boer 2008). In Fig- 
ure 14 we show one such curve, obtained with respect 
to a reference period of 0.5668378 d. As recently re- 
viewed by Sterken (2005), a variable star with a period 
that is linearly changing with time will show a quadratic 
curve in the period-epoch diagram, with the value of the 
quadratic term given by 
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Fig. 15 Average period change rate for RR Lyrae stars in 
Galactic globular clusters. The plusses, stars, crosses, and 
filled gray symbols represent new theoretical results for four 
different metallicity values (as shown in the inset) ; these can 
be very well represented by the analytical formula provided 
in eq. (8) and indicated by a solid line in this plot. 
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where (P) is the average period value during the times- 
pan of the observations, and E is the epoch. A formal 
least-squares fit to the data shown in Figure 14 pro- 
vides a period change rate of ~ —0.1 d/Myr — although 
it should be emphasized that there is no a priori rea- 
son why the period should change linearly with time, 
and therefore no strong reason to believe that the as- 
sumed quadratic law provides an adequate represen- 
tation of the data for RR Lyrae stars (e.g., Rathbun 
& Smith 1997 and references therein). As a matter 
of fact, according to this diagram and the similar one 
presented by Szeidl & Kollath (2000), the rate of pe- 
riod change in RR Lyr has not at all been constant; 
one sees instead an apparent increase in the period at 
JD 2,418,000 or JD 2,420,000, only to be followed by 
a sudden decrease in the period that took place some 
time around JD 2,432,000, with a seemingly oscillatory 
behavior thereafter (interpreted by Szeidl & Kollath 
as being due to the accumulation of random fluctu- 
ations in the period). One way or another, there is 
certainly no strong evidence of a sharply increasing pe- 
riod, as might be expected in the case of a very bright, 
extremely evolved RR Lyrae star that is not strongly 
affected by random period variations. While this by 
no means represents conclusive evidence that RR Lyr 
is not well evolved and therefore overluminous, it is 
at least consistent with the hypothesis that the star 



presents a luminosity that is not dramatically differ- 
ent from that of most other RR Lyrae stars of similar 
metallicity. 



9 Period Change Rates of RR Lyrae Stars: 
The Evolutionary Connection 

As we have just seen in the case of RR Lyrae (Fig. 14), 
more often than not RR Lyrae variables show er- 
ratic rates of period change, presumably due to mixing 
events in the stellar core (Sweigart & Renzini 1979) or 
to the presence of hydromagnetic effects related to the 
conjectured existence of a magnetic cycle similar to the 
Sun's (Stothers 1980). Interestingly, transient magnetic 
fields in the H and He ionization zones have recently 
been suggested to be responsible for the Blazhko effect 
as well (Stothers 2006). Note that the presence of a 
strong magnetic field in RR Lyr (as also required in 
other theoretical scenarios for the Blazhko effect) has 
recently been ruled out by the high-precision spectropo- 
larimetric observations (carried out over an almost 4- 
yr timespan) by Chadid et al. (2004). Stothers (2006) 
speculates, again drawing an analogy with the Sun, that 
these measurements may simply indicate that a strong 
magnetic field may be deeply embedded in the star, be- 
coming rather weak only at its surface (see §3 of his 
paper for extensive references to earlier work on this 
subject). Relatively recent reviews of the Blazhko ef- 
fect in RR Lyrae stars have been provided by Smith 
(1997, 2006) and Kolenbcrg (2002). 

This notwithstanding, given sufficient time coverage 
and large enough samples of stars, one may be able 
to detect, superimposed on the random period changes 
that characterize many RR Lyrae variables, the slow 
period changes that are due to the star's evolution ac- 
cross the CMD. That the period of the star should 
change with time is a very basic prediction of stellar 
evolution and pulsation theory: as a star slowly changes 
its luminosity and temperature, its mean density algo 
changes; according to Rittcr's (1879) period-mean den- 
sity relation of pulsation theory, Py/Jp) « const., the 
period should accordingly change in inverse proportion 
with the density, implying that a star with decreas- 
ing density — which may be due either to a decreasing 
temperature and/or an increasing luminosity, assuming 
the stellar mass is constant — should have an increasing 
period. Since stars evolving to the red and towards 
higher luminosities on the HR diagram are thought to 
be present mainly in globular clusters with blue HB 
types, it follows that positive period change rates that 
can be ascribed to evolutionary effects are expected 
mostly on blue HB clusters. Clusters with predom- 
inantly red or intermediate HBs, on the other hand, 
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should show, on average, period change rates around 
zero, or even negative. 

Some models which clearly show this predicted trend 
of increasing period change rate for bluer HB types have 
been presented by Lee (1991). Recently, Catelan et al. 
(2004) computed a new, extensive grid of synthetic HB 
models, which are also suitable for a reassessment of 
these canonical theoretical predictions. Based on those 
simulations, which wc augmented with an additional set 
of simulations for both extremely blue and extremely 
red HBs, period change rates were computed for the 
metallicity values Z = 0.0005, 0.001, 0.002, and 0.006, 
as shown in Figure 15. Note that these theoretical re- 
sults can be extremely well described by an exceedingly 
simple analytical formula (solid line in Fig. 15), namely 

/AP\ 0.0053 

\aF/ = TT099Z (d/Myr) < (8) 

where C = {B-TZ)/(B+V+TZ). 

The empirical data are also shown in Figure 15, 
overplotted on the theoretical results. Using the re- 
cent Nemec (2004) compilation as a starting point, 
we have added to the sample several additional glob- 
ular clusters, including M2 (Lee & Carney 1999a) and 
NGC 4147 (Stetson et al. 2005), as well as earlier data 
for several globulars as summarized in Table 5.1 of 
Smith (1995) — the latter corresponding to the data- 
points without error bars. These data are listed in Ta- 
ble 4. As can be seen, there is good agreement between 
these evolutionary models and the observations, within 
the empirical error bars. In the case of NGC 4147, 



the timespan of the available observations has presum- 
ably been insufficient to reliably detect the evolutionary 
changes in the pulsation periods. 

10 Conductive Opacities and the He-Core 
Mass at the Helium Flash 

In an earlier review, Catelan, de Freitas Pacheco, & 
Horvath (1996) argued that uncertainties in the con- 
ductive opacities remained that might still lead to sig- 
nificant changes in the computed properties of HB stars, 
especially their luminosities. Since over 10 years have 
passed since the Catelan et al. study, a new look at the 
problem seems especially worthwhile — and this is the 
purpose of the present section. 

To be sure, our approach does not mean that there 
are not other physical processes which are uncertain 
enough as to potentially lead to important changes in 
our understanding of HB evolution; quite the contrary, 
in fact. For instance, the treatments of diffusion, con- 
vection, and mixing, as well as such a key nuclear re- 
action rate as 12 C(a,7) 16 0, remain subject to consid- 
erable uncertainty, and future developments may ac- 
cordingly affect HB evolutionary predictions in quite 
significant ways. Unfortunately, it would not be prac- 
tical for us here to provide an extensive review of all 
of the many different physical processes that play an 
important role in defining the predicted properties of 
HB stars. The interested reader is referred to the re- 
cent reviews by Salaris et al. (2002), Cassisi (2005), and 
Catelan (2007) for critical discussions of several of the 
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Fig. 16 Run of temperature (Ze/Z panel) and density (right panel) as a function of the Lagrangian coordinate mass in 
the core of selected RGB models. Labels next to each temperature or density profile indicate the time (in Gyr) since the 
star arrived at the zero-age main sequence. The top model profile in each panel shows the core profile at the onset of the 
He-flash. The vertical line accordingly shows the mass coordinate where He ignition first takes place. 



more salient uncertainties affecting the computation of 
HB models. 

10.1 Overview 

One of the main ingredients affecting the luminosity 
level of HB stars are the conductive opacities in the 
cores of their progenitor RGB stars. The cores of RGB 
stars are electron-degenerate, so that electron conduc- 
tion is the main form of energy transport. Naturally, 
the more efficient the transport of energy away from 
the He-core, the more difficult it becomes for the core 
to reach high enough temperatures for the onset of he- 
lium burning, with the end result that the He-flash is 
delayed for higher electron conductivities (i.e., smaller 
conductive opacities), the star reaching a higher He- 
core mass and a brighter location on the RGB tip — 
thus also being expected to lose more mass during its 
first ascent of the RGB and thereby settle on a bluer 
post-flash location on the ZAHB. 

The He-core mass M c is also known to have a 
marked, direct impact upon the HB structures. The 
shapes of HB evolutionary tracks are strongly affected 
by the precise value of M c , with a larger M c leading to 
less marked blueward loops on the HB phase (Sweigart 
& Gross 1976) and thereby affecting the predicted frac- 
tion of well-evolved RR Lyrae stars in globular clusters 
with predominanty blue HBs (Catelan 1992). Most im- 
portantly, higher M c values lead to brighter HBs as 



-^gL „ 7.3 - 9.0 (mag/M ) (9) 

(Sweigart et al. 1987; Catelan 1996; Catelan et al. 
1996); therefore, an uncertainty in the value of M c by 
only ±0.01 M Q (or about ±2%) is capable of affecting 
the HB luminosity level by ±0.07 - 0.09 mag. The 
numerical experiments by Sweigart & Gross (1978) in- 
dicate that a change in M c by about 0.01 M Q obtains 
when one reduces the conductive opacities used in their 
calculations by a (uniform) factor of 2. 

The implied level of uncertainty in the HB models 
is important not only for the sake of comparing the 
predicted and empirically calibrated HB luminosity- 
metallicity relationships, but also (and perhaps most 
importantly) in the context of the self-consistent de- 
termination of globular cluster ages from evolutionary 
models of globular cluster stars; of the establishment 
of the RGB tip distance scale; of the determination of 
the helium abundance in globular clusters using Iben's 
(1968) i?-method; and of the placement of astrophysical 
constraints on fundamental (especially particle) physics 
parameters. As to the latter, constraints on the neu- 
trino magnetic moment using bright RGB and HB stars 
were discussed by Raffelt & Weiss (1992) and Castel- 
lani & Degl'Innocenti (1993); those on the axion mass, 
by Raffelt & Dearborn (1987); and those on the num- 
ber of extra dimensions in the Universe by Cassisi et al. 
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Fig. 17 As in Figure 16, but for the run of the Coulomb 
r parameter as a function of mass in the RGB core. The 
vertical line shows the mass coordinate where He ignition 
first takes place. 

(2000). Reviews of the subject and extensive additional 
references have been provided by Raffelt (1996, 2000) 
and Catelan et al. (1996). 

Until the mid-90's, the two main sources of con- 
ductive opacities were the calculations by Hubbard & 
Lampe (1969) and Itoh et al. (1983). To date, most 
authors still use one of these two, as can be seen from 
the recent summary of physical ingredients in different 
evolutionary models compiled by Gallart et al. (2005). 
Yet, as pointed out by Catelan et al. (1996), the the- 
oretical calculations of Itoh et al. in particular are 
not applicable to the physical conditions characteriz- 
ing the interiors of RGB stars, since they are restricted 
to the regime of relatively strong Coulomb intreactions 
(Coulomb parameter T > 2), while in the RGB core 
r < 0.81 — an intermediate regime between the liquid 
and crystal phases considered by Itoh et al. and an 
ideal gas. Since then, however, new conductive opac- 
ity calculations have been presented by Potekhin (1999) 
and Potekhin et al. (1999). 

Figure 16 shows the run of temperature (left panel) 
and density (right panel) in the interior of a metal-poor, 
low-mass star as it climbs up the RGB. These mod- 
els, which refer to a 0.8 M star with [Fc/H] = -2.3, 
[a/Fe] = 0.3, and Y = 0.2484, were kindly provided by 
D. A. VandcnBcrg (2005, priv. comm.). Nine "snap- 
shots" are taken over an interval covering the final 
56 x 10 6 yr of the star's evolution in the RGB phase; 
as can be seen from this plot, during this time interval 
the mass of the He core increases from about 0.28 M Q 



all the way up to its final value, of order 0.5 M . Note 
that the He-flash actually takes place at a mass value 
just below 0.26 M Q , due to the fact that the density- 
dependent neutrino energy losses lead to an efficient 
cooling of the innermost regions and therefore to a tem- 
perature inversion in the RGB core, as first shown by 
Thomas (1967). For the sake of simplicity, in what fol- 
lows we shall assume a pure helium plasma. 

Figure 17 shows the variations in the Coulomb T pa- 
rameter [see, e.g., eq. (21) in Catelan et al. 1996], which 
measures the strength of the electrostatic interactions 
in the plasma, throughout the He core as a function 
of the star's evolutionary status. Clearly, this diagram 
fully supports the assertion that T < 0.81 in the RGB 
interior. This accordingly confirms the caveat that one 
should avoid using the Itoh et al. (1983) results in evo- 
lutionary calculations of low-mass stars. 

Since the more recent Potekhin (1999) and Potekhin 
et al. (1999) calculations have superseded the Itoh et 
al. (1983) results, the preceding argument leads natu- 
rally to the following question: are the Potekhin results 
themselves fully applicable to the conditions character- 
izing RGB interiors? 

In fact, not quite: Potekhin et al. (1999) explain 
that a fit to the static structure factor that is valid 
over the range 1 < T < 225 was employed in their 
calculations. Therefore, the possibility that a problem 
similar to the one affecting the Itoh et al. (1983) results 
is present also in the case of Potekhin's calculations 
cannot be excluded; a fully self-consistent calculation 
for 0.1 < r < 0.8, which is the relevant range for the 
typical RGB interior (Fig. 17), would prove of great 
interest in the context of precision models of RGB stars 
and their HB progenitors. Potekhin et al. do note 
that Coulomb logarithms in the transition domain from 
weak (r <C 1) to strong (T > 1) ion coupling could be 
calculated using the Boerker, Rogers, & DeWitt (1982) 
formalism, but unfortunately such a formalism was not 
applied in their calculations. In fact, it appears that 
the reason why such a range of T values has not been 
given great emphasis in such calculations as those by 
Itoh et al. or Potekhin et al. is that such studies are 
generally carried out primarily with applications to the 
more extreme conditions characterizing neutron stars 
and white dwarf cores in mind, the interiors of RGB 
stars having received comparatively little attention. 

This, in fact, also brings to mind the question 
whether other approximations to physical situations 
that are rather more extreme in neutron stars and white 
dwarf cores may not have been adopted in these latest 
calculations that could render uncertain their applica- 
tion to the conditions characterizing RGB cores. In 
this sense, the most obvious ingredient is the electron 
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Fig. 18 As in Figure 16, but for the ratio between actual 
temperature and the Fermi temperature in the RGB core. 

degeneracy level. Figure 18 shows the run of the ratio 
between the actual temperature and the Fermi tem- 
perature for the same models shown in the previous 
plots. As well known, for T <C TFermi, strong degener- 
acy is present, whereas T 3> 7F orm ; implies a classical 
regime where degeneracy effects are unimportant. This 
plot indicates that RGB interiors fall in an intermedi- 
ate regime between strongly degenerate and classical. 
In other words, partial degeneracy better describes the 
status of the RGB interior. Note, in particular, that at 
the mass coordinate where the He-flash takes place, the 
temperature ranges from about 10% to 45% the Fermi 
temperature. Clearly, the effects of partial degeneracy 
should be fully taken into account for conductive opaci- 
ties that are applicable to the conditions characterizing 
RGB cores to obtain. How does this compare with the 
Potekhin (1999) and Potekhin et al. (1999) analyses? 

As clearly stated by Potekhin et al. (1999), their 
analysis "is limited by the condition T <C Tpermi" — 
in other words, strong degeneracy is assumed in their 
calculations. Unfortunately, this cannot be expected to 
provide particularly reliable results for the RGB inte- 

97 

nor. 

The above discussion refers to the contribution of 
electron- ion (ei) interactions to the conductive opac- 
ity. An apparently more subtle effect is related to the 



27 On the other hand, Potekhin (2005, priv. comm.) points out 
that the technique of "Fcrmi-Dirac averaging" has been adopted 
in the region of partial degeneracy, which at least improves some- 
what the situation compared to what would obtain from straight 
application of the strongly degenerate results to partially degen- 
erate conditions (see Potekhin 1999 for more details). 



Fig. 19 As in Figure 16, but for the ratio between the com- 
ponent of the conductive opacity due to electron-electron 
interactions and the total conductive opacity, according 
to the calculations by Potekhin (1999) and Potekhin et 
al. (1999), in the aforementioned RGB core models. The 
electron-electron contribution is clearly a major component 
of the conductive opacity for the conditions characterizing 
the cores of RGB stars. 

electron- electron (ee) interactions. Since these are usu- 
ally of little importance in the more extreme condi- 
tions characterizing neutron stars and white dwarfs, the 
Potekhin (1999) and Potekhin et al. (1999) calculations 
have not treated this component in as much detail as 
in the case of the ei interactions. The ee component 
was indeed computed assuming complete degeneracy, 
with no corrections whatsoever due to partial degen- 
eracy effects — and are accordingly expected to be ac- 
curate only to within a factor of ~ 2 (Potekhin 2005, 
priv. comm.). However, it must be noted that ee effects 
are a major component of the conductive opacity in a 
He plasma in the conditions characterizing the RGB in- 
terior (Hubbard & Lampe 1969; Catelan et al. 1996). 
This is shown in Figure 19, where the ratio between 
the component of the conductive opacity due to ee in- 
teractions and the total conductive opacity, according 
to the Potekhin calculations, is displayed. One sees that 
the ee contribution reaches a minimum of about 28% 
of the total conductive opacity close to the mass coor- 
dinate where the He-flash takes place. The ee contri- 
bution increases further as the He-flash is approached. 
Evidently, an error by a factor of two in the ee contri- 
bution to the conductive opacity implies that the total 
opacity is uncertain by a comparable factor. Accord- 
ingly, greater attention should be devoted to this crucial 
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Fig. 20 Electric conductivities: comparison between lab- 
oratory experiments (black symbols with error bars, from 
Stygar et al. 2002) and the expected theoretical values ac- 
cording to Potekhin (1999) and Potekhin et al. (1999) (gray 
symbols). 

physical ingredient in future calculations of conductive 
opacities for RGB stars. 

10.2 Laboratory Experiments 

It is obviously of great interest to check how the 
Potekhin (1999) and Potekhin et al. (1999) opacities 
compare with experiments conducted in the laboratory. 
Such a comparison, in the case of electrical conduction, 
has recently been carried out by Stygar, Gerdin, & Fehl 
(2002), using laboratory data from Mintsev, Fortov, & 
Griaznov (1980), Benage, Shanahan, & Murillo (1999), 
and Benage (2000). Their results are shown in Fig- 
ure 20, where the electrical conductivity is plotted as 
a function of the Coulomb T parameter. These com- 
parisons clearly reveal a very complex pattern, with 
the models systematically overestimating the conduc- 
tivity (i.e., underestimating the conductive opacity) for 
r > 1.9; presenting good agreement with the empir- 
ical results for 1.3 < T < 1.7; and again overesti- 
mating the conductivity for 0.85 < T < 1.25. It is 
unclear at present what the reasons for this behavior 
are, though Potekhin (2005, priv. comm.) points out 
that the laboratory plasmas, when heavy elements are 
present (the above experiments were carried out using 
xenon and aluminum), are almost always partially ion- 
ized, and it is very difficult to determine, and therefore 
account for, the actual degree of ionization, because of 
the plasma nonideality. Accordingly, more empirical 
tests of the Potekhin results using low-Z gases, espe- 
cially He, would be highly desirable. 
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Fig. 21 Relative difference between the conductive opacity 
computed by Itoh et al. (1983), with quantum corrections by 
Mitake et al. (1984), and by Potekhin (1999) and Potekhin 
ct al. (1999) in the aforementioned RGB core models. The 
inset shows a blowup of the region indicated by a rectangle 
in the upper left of the plot. 

10.3 Comparison between Different Prescriptions 

It is instructive to compare the Potekhin (1999) and 
Potekhin et al. (1999) results with those from the Itoh 
et al. (1983) and Hubbard & Lampe (1969) analyses. 
This is done in Figures 21 and 22, respectively, which 
show the relative differences between the conductive 
opacities computed by Itoh et al. and by Hubbard 
& Lampe with respect to the Potekhin results for the 
physical conditions characterizing the RGB interior cal- 
culations. In the Itoh et al. case, we have included 
the quantum corrections by Mitake, Ichimaru, & Itoh 
(1984), which however are rather small over the region 
of interest. The agreement between the Itoh-Mitake 
prescriptions and Potekhin's is quite poor, with the for- 
mer systematically overestimating the conductive opac- 
ity except very close to the core, where the conduc- 
tive opacity is instead underestimated by about 25%. 28 
On the other hand, the Hubbard & Lampe calculations 
reveal a better agreement with the Potekhin results, 
the maximum overestimate of the conductive opacity 
by Hubbard & Lampe compared to Potekhin not ex- 
ceeding about 55%. In the central core, again the older 



28 According to Potekhin (2005, priv. comm.), this discrepancy 
can largely be ascribed to the fact that the Itoh et al. (1983) 
calculations are limited not only to strong coupling, but also to 
strong electron degeneracy — compare with the discussion in the 
previous footnote. 
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Fig. 22 Relative difference between the conductive opacity 
computed by Hubbard & Lampe (1969) and by Potekhin 
(1999) and Potekhin et al. (1999) in the aforementioned 
RGB core models. 

calculations underestimate the conductive opacity com- 
pared to the more recent results, albeit by a maximum 
of less than 20%. We conclude that the Potekhin re- 
sults are in much better agreement with the canonical 
prescriptions by Hubbard & Lampe than they are with 
the results by Itoh et al.. 

While updated conductive opacities that are entirely 
suitable for the conditions characterizing RGB interiors 
are not available, it is clear that one should avoid using 
the Itoh et al. (1983) results in calculations of low-mass 
stellar evolution. 

10.4 Latest Developments 

After the first version of this paper had been completed, 
new conductive opacities were published by Cassisi et 
al. (2007) — in essence, a revised and updated install- 
ment of the Potekhin (1999) opacities. These calcula- 
tions, while still not perfect, already do take several of 
the shortcomings noted in the previous sections into due 
account, and should accordingly be strongly preferred 
in state-of-the-art evolutionary computations. 

§2.4, and especially Figures 3 and 4, in Cassisi et al. 
provide the results of a comparison between the new 
conductive opacities computed by Cassisi et al. and 
the earlier studies by Hubbard & Lampe (1969), Itoh 
et al. (1983), and Potekhin (1999). Important differ- 
ences are found with respect to the results from all 
previous studies, especially those from the Itoh et al. 
team. According to Cassisi et al., the differences are 



mainly due to the following: i) Deficiencies in the ear- 
lier treatment of strongly coupled and relativistic plas- 
mas (Hubbard & Lampe); ii) An inadequate extension 
towards the T > Tp regime (Itoh et al.); iii) A neglect 
of ee scattering (Potekhin). 

When the new conductive opacities are used in evo- 
lutionary computations, Cassisi et al. find that the 
changes, with respect to the results based on the ear- 
lier Potekhin (1999) prescriptions, to be relatively mod- 
est. More specifically, they find: i) A reduction in the 
core mass at the He-fiash by about 0.006 M Q , irrespec- 
tive of metallicity; ii) An increase in My (ZAHB) at the 
RR Lyrae level, by an amount ranging from 0.06 mag 
at Z = 0.0001 down to 0.04 mag at solar metallicity; 
iii) An increase in the HB lifetime, by about 5-6%. 
While relatively small, these systematic effects must 
obviously be properly taken into account in future pre- 
cision studies involving HB stars. 

11 Epilogue 

The study of horizontal branch stars encompasses 
and/or has direct implications for many different 
branches of astrophysics. Without first looking into 
several of these areas, one might rush into the con- 
clusion that the study of HB stars is a frustratingly 
limited affair. Nothing could be farther from the truth. 
As we have seen, in order to properly understand HB 
stars and appreciate their astrophysical implications, 
one must dwell on the physical processes (both canon- 
ical and non-canonical) that control the helium core 
flash in RGB stars; the physical processes that lead to 
mass loss in RGB stars; whether (and how) deep mixing 
takes place in RGB stars, both before and during the 
helium flash; how angular momentum evolves with time 
in low-mass stars; where, and how fast, radiative levi- 
tation and diffusion effects in moderately high-gravity 
stars "kick in" ; how dramatically non-solar abundance 
ratios (again as a result of radiative levitation and grav- 
itational settling) may affect model atmospheres and 
the photon output as a function of wavelength; how 
both radial and non-radial (both p and g) pulsation 
modes can be excited, and then evolve with time; how 
theory and observations of radial pulsators may place 
constraints on the physical parameters of stars, and 
thereby help determine the formation history of galax- 
ies and the extragalactic distance scale; how astero- 
seismology of non-radial pulsators may help reveal the 
interior structure and thereby constrain the evolution- 
ary history of stars; how the ultraviolet flux from giant 
elliptical galaxies and bulges of spirals come to being; 
how the stellar populations of resolved and unresolved 
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galaxies evolve with time; and so on and so forth. In 
summary, the study of HB stars is a challenging and 
far-reaching intellectual adventure, and we hope that 
the present review will have helped unveil some of the 
reasons why this is so. 
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